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This thesis has dealt with various approaches to modelling 
2-D data fields using state-space formulations. Computer 
Simulation of these models has been carried out to generate 
Simulated 2-D data which could then be used for various 
other signal processing operations. An interesting develop- 
ment that has resulted from this study is that of adaptation 
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Lo. LNDTRODUCTION 


ATHE MAIN IDEA 

Image processing by nonoptical means has received exten- 
Sive attention in the last few years. Several books and many 
papers have been published that have established nonoptical 
image processing as a viable area of research. A large 
portion of this research emphasizes the linear processing 
of images for two main reasons: 1) Many image processing 
tasks are linear in nature. These tasks include image enhance- 
ment, image restoration, picture coding, linear pattern 
recognition, and TV bandwidth reduction. 2) There are many 
known linear techniques that may be brought to bear in the 
treatment of linear image processing. These techniques include 
transform theory, matrix theory, filtering, signal modeling, 
etc. Several common operations involved in image processing 
include transfer function concepts, partial difference (recur- 
Sive) equations, and convolution summations. For example, 
Vander Lugt [Refs. 1,2] has presented an extensive development 
of linear optics based on transfer functions. The transfer 
functions relate the two-dimensional Fourier transform of an 
output image to that of the input image. Complex optical 
systems are easily described by combinations of transfer 
functions that correspond to individual components of the 


Optical system. 


Partial difference equations are used by Habibi [Ref. 3] 
to describe a model for estimating images corrupted by noise. 
The model corresponds to a two-dimensional extension of Kalman 
filters.  Convolution summations are discussed by Fryer and 
Richmond [Ref. 4] in work that involves simplifying a two- 
dimensional filter to a single dimensional filter. 

The time-discrete state-space model offers great 
utility in the formulation and analysis of linear sys- 
tems. Linear systems that are described by transfer functions, 
difference equations, or convolution summations are formulated 
into a state-Space representation. Once formulated, many 
known techniques may be applied to systematically analyze 
the model. Consequently, the state space model is a general 
and powerful tool that is used to unify the resedi and the 
study of time-discrete linear systems. 

This thesis develops the discrete model of Roesser [Ref. 

5] for linear image processing which closely parallels 

the well-known state space model for time-discrete systems. 
Because it is parallel, many of the concepts that are known 
for the temporal model may be carried over to the spatial 
model. This is done by generalizing from a single coordinate 
in time to two coordinates in space. The spatial model will 
hopefully have some of the same utility for the study of two- 
dimensional linear systems as the temporal model for one- 
dimensional linear systems [Ref. 3]. However, not all of the 
properties of one-dimensional systems carry over into the multi- 


dimensional case. 


One of the fundamental problems involved with recursive 2- 
Dimensional systems is that the order of the System (recursive 
memory) is not the same as the number of initial conditions 
(boundary conditions). In one-dimensional systems these are the 
same. Temporal systems are inherently nonanticipatory and are 
peuEEreated as such for the sake of physical realizability 
in real time; whereas spatial systems do not have causality which 
15 an inherent limitation. That is, an image processor may have 
right to left dependency as well as left to right dependency. 
Finally it is noted that stability criteria in one-dimensional 
recursive systems become much more difficult when carried over 
to the multidimensional case. 

Causality is built into the temporal state-space model if 
an initial state is assumed to be fully specified. In order to 
establish a close parallel for the spatial model, the same 
built-in causality will be intentionally assumed, despite the 
fact that causality is not necessary for physical realizability 
in real space. Such an image processor is said to be unilateral. 
If the constraint of causality is removed, then the image proces- 
Sor is said to be bilateral [Ref. 5]. Concepts that are 


developed in this thesis for the latter case are: 


1) Formulation of the state space model of Roesser. [Ref. 5] 
2) The definition of state transition matrix. 

3) A resulting computer program based on the above model. 

4) An investigation of the class of 2-Dimensional transfer 


functions defined by this model. 


5) Derivation of a general response formula. 


6) Extension of Roesser's model of state variable equations 
ἕο encompass a largemyclass Of fran EEr COC 


7) Adaptation of the 1-D "SSPACK" program to produce 2-D data. 


I STATE SPACE REPRESENTATION 

Toward the end of the 1950s, the concept of representing 
a discrete system by a set of first-order difference equations 
became a standard tool of the research engineer. These tech- 
niques have since become generally known as state-space repre- 
sentations. Such representations have become increasingly 
important during the intervening yearns because the .:-.; 
one to carry out a meaningful system design entirely in the 


discrete-time domain (in comparison to popular Z-transform 


methods). That this is important follows basically from these 
incor: 
l. The system may be nonlinear so that transformation 


methods are not directly applicable. 

2. Time-domain concepts often give one a better insight 
into the analysis and synthesis of the system (fre- 
quently with the aid of a digital computer). 


3. Cases in which the initial conditions are non-zero may 
be handled straightforwardly. 


A state space representation of a system differs from the 
conventional representation. In a conventional representation 
only the relationships between the input and output signals need 
be known. On the other hand, the state-space representation 
gives a total description of both the internal as well as the 


external signals of a system. 
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C. STATE-VARIABLE REALIZATIONS--THE CONCEPT OF STATE -- 

ο oO limeame systems theory and@eontrol theory, the 
concept of a filter state has played an important role. 
Basically the filter state at any point in time contains 
all the information necessary to compute the remainder of the 
πι output signal, given the “imput signal. Ome dimensional 
ου ος Ἱπρος, single-output filter realizations based on a 


state variable model can be written in the forn: 


x(k+1) = Ax(k) + Bu(k) ει) 


y(k) = Cx(k) + Du(k) (1.19) 


G rorm relates the input u(k) and the output y(k) through 

a state vector x(k). “he state vector evolves in time accord- 
mee omequatlen (1.14). The matrices A, B, and C and 1x1 matrix 
Ευ ern the exact form of the input-output relationship. 

(In general these matrices may vary with the index (k) and 

the input and output signals may be vectors as well.) Quite 
often the components of the state vector are taken to be the 
constants of the pL delay operators in a flowgrapn represen- 
EIN Of the 1=D filter. 

A classic problem in state-variable theory representation is to 
find the matrices A, B,Cand D which willrealize a particular system 
function H(z) with a minimum number of state variables. A 
Similar approach may be taken to develop a 2-D state-variable 


model. 
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A 2-D discrete system may be defined as a mathematical 
abstraction which utilizes three types of variables to repre- 
sent or model the dynamics of a discrete-time process. The 
three variables are called the input, the output, and the 
state variable. The input variables u(i,j), serve as external 
forces which influence the dynamics or motion of the system. 
The output variables y(i,j) are the characteristic variables 
which are directly observable (measurable) by the external 
observer. The state variables x(i,j) characterize the internal 
dynamics of the system and are to be selected according to the 
following rule. 

These variables are formulated in such a manner that, if 
one knows the values of the present state variables x(i,j) along 
with the values of the input variables u(i,j) then the output 
variables y(i,j) and the next state variables x(i,j) are com- 
pletely determined. Moreover, the number of state variables 
used in a state-space representation must be minimized. A ` 
state-space representation may be visualized in block diagram 


form, as shown below. 


ie E γη} 
Si AT ο 


X). , XG3)...X G) 


yid) 





Figure 1.1 


2 


In Figure l.l1,m-inputs, p-outputs and n-state variables are repre- 
sented. However, we willbe mainly interested inthose systems which 
have one input (m = 1) and one output (p = 1). It is important 
to note that the input and output variables appear external to 
the system, while the state variables are generally internal. 

The different input variables will be represented by the 


input vector u(i,j) where, 


ΠΡ) 


the output vector y(1,j]) where, 





ana the state vector x(1,]) where, 





For a given process the state space representation is not unique. 
However all such representations have one characteristic in 
common for a given system, namely the number of elements n is 


K = y to as the order of the system. 


ES 


II. ROESSER نت‎ 2110 ΠΟΘΕΝ 


A. THE FRAMEWORK 
An image is a generalization of a temporal signal, in that 

it is defined over two spatial dimensions instead of a single 
temporal dimension. Consequently, two space coordinates i and 
j take the place of time, t. Also, two-state sets are intro- 
duced to replace the single-state set. The following defini- 
tions are made by the model: 

1 An integer-valued vertical coordinate; 

j An integer-valued horizontal coordinate; 


{R} A set of nj, real vectors which convey information 
1011312011311 s 


{9} A set of n, real vectors which convey information 
vertically; 


{u} A set of m real vectors that act as inputs; 
(v) A set of p real vectors that act as outputs. 


A Specific image processor is then defined as 6-tuple 


)٣٩(, ۱5 TUE 2 τα E 


where f is the next sto cm Une E 


foad iR απ ο. ο e 


and y is the output "σπα 177 


q: 11 ٤١ ΠΗ ΤΙ ΟΡ ΘΙ 
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Now since f and g are to be linear functions, they may be 


represented by the following matrix equations: 


R(itl,j) = A,R(i,j) + A,S(i,j) * B,u(i,j) 
s is jp a= AjR(i,j) + A,S(i1i,j) + B.u(1,3j) (IPIE S 
NE Η; CIE Ec sg SEEMS O 


Άη, Ay, Az, Ayr Bi, Boy Cir Cos D are matrices of appropriate 
dimensions. Boundary conditions R(0,j) and S(i,0) and also the 
input u(i,j) are externally specified. In the next section a 


computational rule is obtained that uniquely OSO ses the 
ERE 370 51 ,(( and also the output y(i,j) (formi,j > 0) 
from the boundary conditions (such as all zero). The equations 
pasduseca set of output vectors from the input vectors. 

This formulation is general so that any discrete linear 
image process may be so represented. Notation is condensed 


somewhat by introducing the following matrices and vectors: 


R8 - ee 1 C = [C C] 
ο By 
ru) R(i*l,j) 
ο = TE = 
S(i,j) S(i,j+1) 
η"; - APO cL Bul.) 
pO) αρα πλ δι ντ) 


ΠΠ ΕΕΞΡΌΝΒΕ FORMULA 


A state-transition matrix A is defined as follows: 


15 





MÀ 1 η. 
Then exponentiation A 'J is defined as, 


alriat Jnl Ca oon‏ " ,لحد لال 


= a’ J = ο for ۱۹٣۳۷ eee 


Examination of this definition bears out that it is an effective 
recurSive definition of ΠΠ for integer values of 1 and j such that 
either i > 0 or j ο ο ος IT Es OE دن‎ sns the defini- 
tion of the time-discrete state-transition matrix As = A ee 

It now remains to be shown that this state transition matrix 
E may be used in expressions for the response of the model 
in terms of the inputs and boundary conditions. The term 
boundary conditions is used here to refer to the states along 


the edges of the model. Specifically, the set of boundary condi- 


tions consist of R(05;j 6 BR n πο οι ο ο ο Ὁ- 


C. CHARACTERISTIC FUNCTION OF A MATRIX 
If the primary inputs and outputs are dropped in the model 
equations (II.1), a representation arises for the state behavior 


of the system having the form 


ΤΠ 


A,R(i,j) = A,S(i,j) 
(TI : 9 
S 1 ٢ 


II 


A4R(i,j) + A,S(1,J) 


These equations are useful in the development of a form for a 


two-dimensional characteristic matrix of A. Operators are 
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first introduced that advance a particular coordinate of 
their operand. 

Definition: Let E be an operator that has the effect of ad- 
vancing the vertical coordinate or the first subscipt of the 
function upon which it is operating. Likewise, let F be an opera- 
tor that has the effect of advancing the horizontal coordinate or 
EcesDdosubscript of the function upon which it is operating. 


The effect of these operators on the state vectors is: 


R(it*1,i) 


ERT 


S(i,j+l) = FS(i,j) 


The state equations can be rewritten using these advance 


operators. 


I 
C 


ΤΡ - A48(i,]j) 


II 
C 


-AjR(,j) τ (FI-A,)S(1,3) 


These equations are equivalently represented in the overall 


matrix form. 


(EI-A, ) _ A 
l : T 


-A (PTER 


3 4? 


The above equation represents a system of homogeneous equations in 
the elements of T(i,j). If the system is to have a non-trivial 
solution for T(i,j) then the transformation represented by the 
matrix must be singular. The above matrix is said to be the two- 


dimensonal characteristic matrixof the partitioned matrix A, where 
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The characteristic matrix of A is denoted cm(A) and may be 


represented as 


cmA πε prt? + pporl - A 
where, 
j 
TONO ο Í! O 
,ل1‎ E τομ and 0,1 2 e 
REN 0 I 1 


Now since cm(A) must be Singular, 1ts determinant must ος ο πα 
to Zero. | cm (A) | = 0. If E and F are replaced in the above by 
general indeterminates x and y respectively, the result is an 
expression called the two-dimensional characteristic equation 
for A. The determinant of cm(A), and x and y replacing E and 

F. is called the two-dimensional characteristic function of Ne 


matrix and is denoted by 
πο EU MM 


f(x,y) will be a monic polynomial in x and y with degree n4 


in x, and degree n, in y, Where n, is the dimension of R and 


Jà 
n» is the dimension of S. f(x,y) has the form 
κ ο... ) ) 2 yey? 
(0,0)<(i,j)<(n n>) 1 
where a. . denotes elements of A anda = 1. 
وول‎ Τη «Πο 
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D. CIRCUIT ELEMENTS AND THEIR REALIZATION 


EUS Consider tie single 2-D IIR filter transfer function 


glven by: 
8 -1 E -2 -1 
٢ b z 
H(z;,2,) = ΠΡ 01-2 | nee TT 5 0 ΒίΖι 2») 
E = E τ -1,-1 TANE: Ez E) 
1 21071 730172 731174 75 82121 72 و د‎ 
A simple block diagram for Η (71/22) follows. 
y (ii) 






W (id) 





1316 2-7 


Mepa Ee Signal u(i,jJ) flows through a filter corresponding 


peemumerdtor transfer function B(z,,242). The resulting 


1 
Signal 1s added to the signal-w(i,j) to produce the output 


) 


signal y(i,j). The denominator transfer function 1-A(2,;2, 


is realized by the feedback loop containing A(21;25). 
Since we are dealing with two dimensions, there are two 
fundamental shift operators which may occur along a signal flow 


path, the horizontal shift operator indicated by zl and the 
vertical shift indicated by 24 [we shall omit from consideration 
the inverse shift operators Z4 and zl. 18 1٢۵ اه کت‎ Cases -of prac- 
tical interest they can be eliminated by multiplying both the 
numerator and denominator polynomials of H(z4,25) by the appro- 
priate powers of η, and 2 Let us look at a signal flowgraph. 


representing the numerator polynomial: 


go 


= b + b m = É -1 -1 -2 -1 


B(Z, 2.) 00 1021 T bo 122 BD s E 


i 


which is shown in Figure 2-2 below. 





Figure 2.2 


Note the chain of two es Operators descending on the left and 


the single σα, Operator ascending on the right. The nodes 


along these two vertical paths are connected by branches with 
the appropriate gains. If we label tne nodes in both وو‎ 
Chains and the F chain 0,1,2 and so on, fromthe tan down, 


i 


the ith node in the 25 chain is connected to the jth node 


J 
Similarly the signal flowgraph for the polynomial A(z4,2.) 


in the zm chain by a branch withwascammeractor Of b, 


is shown in Figure T:mi 


UN eut 


OUTPUT 





Figure 259 
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Since there is no Ang term, there is no direct connection be- 
tween the input and output nodes of this Signal flowgraph. Thus 
any path from the input node to the output node will encounter 
at least one Ta Or ΖΝ shift operator. 

At this point it iS appropriate to discuss realizations for 


the two shift operators zl and T At their simplest level, 
the shift operators merely select the "previous" S-tuple value 
in the horizontal or vertical direction. When the input to a 


Ze Seerator is the S-tuple u(i,j) the output will be R(1-1,3). 


0 1115371 7 for a Z, Operator "ΕΠ ΙΙ ΤΕ usi hen 

1 تت‎ ορ. τε Ria, 7) Or S(1,J). Consequently a realization of 

either shift operator must embody the appropriate amount of memory 

to retain the "previous" S-tuple in the appropriate direction. 
Interestingly enough, in the more general case where the numer- 

ator and denominator polynomials are considered jointly, the state 


variable realizations based on conventional signal flowgraphs may 


ου TS minimal im the sense that the transfer furnction can be 


realized with fewer coefficients. Consider, 
ES =. وات‎ 
ποιοι» 
Η(Ζη Z.) = err ET CT GNI OE 
η. JS NEST 2 


The corresponding signal flow representation is shown in Figure 


Ze Selow: 
GM 





Pigure 2.4 
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E.- ANALYSIS OF ROESSERGSSMODEI 


Recalling from page 14 the equations of the model are: 


٥ M E A,R(1,)) + A,S(1,J) + B,u(i,)) 
۹۱1 E A,R(1,j) + A S (1,3) + Bou(i,j) 
y(i) = C R(I J) πι ο ο τῇ -υις 


al D 





ΔΙ, A^, .په‎ Agr Bi; B5, Ci: C5! D are scalars or matrices of 
appropriate dimensions. 
R(it*l,j) A, Εν η) Bi 
= + 1171 (le) 
Sur. t A, ene.) B. 
RJ) 
xu gu Em [e QM : ΠΤ (11.9) 
T 2 ΜΉ 
ΕΞ 
κι τη; ᾿Ξ AJR(i,j) + Δ25(1, 1) + B,u(i,j) 
SEL J) ΙΡ AjR(i,j) + A,S(1,J) + Bou(i,}) 


And taking Z transforms: 


2182122) E A, R(2, 725) + 52952122) + B,u(2,,2,) 
22357122) Ξ A,R(zZ, +25) + A,S(Z125) + Β.α (5122) 
Ε(21 22) 
Y (21:25) - [σι 623 + Ώα (2122) )11 G 
S(z4:25) 


2 


ος 


2182122) = Άη Ἀ(Ζη 2ο) - A,S(2 12,5) = B,u(Z,/2,) 
Z5S(21,2Z,) = A R(Z} 2ο) - Ά 5 (71,72) = B,u(Z,,2,) 
er 
R(Z, 725) [2, -A,] = Α235 (7122) = B,u(z,,2,) 
R(Z,/25) [-A3] - [Z,-A]S(z,,2,) = B,u(Z,/25] 





Or I 
ROZ: Za) B 
ΙΟ Z 1 UCM 
S (21/25) Bo 
er 
where Z; = Z, I and Z 
R(Z} 72) 





T Du (24 20) 


Or 





y (21725) 
l (2 525) 
1 2 u(Z,+Z,) 
-1 
0 A, 1 As 5, 
- iC, Co] - -- = کاک‎ = 
2» η 5, 
eet E 2 


The submatrix 21 1S Simply z, times an identity matrix of the 


1. 


appropriate size. Similarly 22 15 Z times an. identity matrix. 


The objective of the state variable realization procedure is 


to find the matrices A, B, C, and D which yields an F(Z; ,Z ) 


2 


that equals or approximates a desired system function H(z_,z 


1 2) 
In essence, the eguations of Roesser represent an implementation 
for which a design algorithm must be found. One choice for 
the state variables is the output signals from the shift 
operators. 

Thus R(i,j) is a vector containing the output signals 
from the 21. operators and S(i,j) contains the output signals 
from the 225 Operators. (Note that the output signal of a 
shift operator signal path is not necessarily the same as the 
nodal signal at the node to which the signal path points.) 
If a state variable corresponds to the output of a shift operator, 
the next value of that state variable must correspond to the 
input of the shift operator. To obtain the submatrices Άι, 
Ay, A^, Ay in equations of Roesser, we write the input signal 


of each shift operator in terms of the outputs of all the 


24 


shift operators, taking care to include all shift-free paths 


from oupuwt to input (seethe following flowgraph). 








m s NAI] 
U (i) Mu) 
Figure 2.5 
Expanding the form of Equation II-7, page 24, yields: 

| = 
ېت‎ | 0 . B, 
Γρ ο = [C x por. = πρ...» 
٣ 2 Ji 2 0 7 5 
mp 2 
ex 1 
A qat A adj A 


2 


p c] 1 













1 2 | (2 A DOS WP = 
E 
πα 7 )B 
A cqui P 
απο a) 7; 
= [C4 
A,Bi * (217A1) Bj 
(2) -Aq) (25784) -A,A, 
Or 
brin s C4 (257A4) B4 + C1A5B5 + CAB} + C,(2,-A,)B, 
= er ae. جيسيم‎ 
٢ (24 7A]) 125 7A4) -A4À4 
Or 
πα με} د د هوي دل هل اور ددد‎ 
, = E a 
L2 2122 7542] 7 A155 7 A4, FAA, 


(C; 4B, *C)ÀA,B, -C,B)A,;-C Bu A4) c (CB524 *C, B, 27) 
- A Z 


(A, A47 A5À4] πι > Se 


(Lia ey 


Equating equation {11.5} ντε τ. ου ce N yields 
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(C4 A,B,*C,A,B, -C,B, A, -C B, A, ) 4C B, Z. 4C, B, Z., 
| = - 
(AjA,7A,A.) - AzZ) - AZ) + په‎ 


i | N 


Ti si يه‎ 
91921 ΤΌρΙ22 تو لل‎ 
l-a z` l-a 5 و‎ 1-1 

ο το E 


For this example, Z1 = Z); 2» - Z5, all of the coefficients on 
the left hand side are scalars. Equation terms of equal powers 


eee. and z 


1 2' 
0125216284 512281” 61818, - Dia = 0 ` πε ο] 
CoB = bio (II.10) 
CB = bg ΠΠ 
ALA, = AAA. Ξ 3 ΜΕ 1ο 
A, = ap ILS 
AJ = aol (II.14) 
ai] 7 =l ) G) 
From these equations, assuming that B4 - B, - lo τς ποιο ο that: 
s Πο 
x NT 
ΠΣ iom 
DE 


D 


From Equation ph h> 


A.A, - A.A = l = =a 


n Daas ie 
سک لا‎ ο απ η ο οι) 
A.A = a E 


ie. SOS 


Let A, and A, take on particular values p and q respectively, 


or 
pa = ai amol EC 
From Equation (II-6): 


CAB; + CA ο η a 


ὃς, 


b p + b - b - 591810 = D11 = 0 (II.DEB 


01 108 10301 


Substituting Equation (IL. lG P i oren pho ΠΙΟ 


3 22 Ge 
E 1T 2000) m 


01 d 108 πω - ο 


-b τ οι ΠΠ 


10501 


ΟΥ 


2 8 


2 
i O LIO“ * (bgja8,,*bg1a84,989,) = O0 . (II.18) 


The results are just the same as in [Ref. 8]. After the com- 
parison between Roesser's model and the 2-D IIR filter, 


described by Equation (II.3), we have: 


f ao 
11 ο το t οχι - 21529] 
acer Y! bio - (b, 9294 *b94 34 9*b, 1) q* (ba; هم‎ | 
U 01 
| og esq 
SM - b πι πο) 
rom» 
By = 1 
B, - 1 
D = 0 


The foregoing equations relate the coefficients of the 2-D 
transfer function to the terms of the system matrices of the 
Roesser model, Equation (II.1). 

Kung et al. [Ref. 2] have shown that the following state 
variable equations, which use only two shift operators, will 


also realize H(z1:25). For the foregoing example, 
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R(itl, 7) a 


| ° 
8 ROI ile 
= ELE + ú Om) 
5(4, 131) q ‘aig [|5{1,}) 1 
Lu 
G P D 
D. S (i,j) 
OL 
= =] 
zd I 1 
0 Lu 
H(zj,25,) = χο 2573 E a : 
d 2 o 


We can construct a signal flowgraph with only two shift 


Operators. It is an equivalent figure to that on page 25. 


Kung et al. [Ref. 2] have also shown that state-variable 


realizations of the form of the equations above may be 
generalized for any system function H(ZJ Z.) which satisfies 


the following three conditions: 





Xos 
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1) The constant term in the numerator, boo = 0, must be 
zero. 


2) The largest powers of z in the numerator and 
denominator polynomials, must be equal, and 


3) The largest powers of 2 in the numerator and 
denominator polynomials must be equal. 


There is one potential difficulty with state variable reali- 
zations of this type. The nonlinear equations defining p and 
q may result in complex values for these constants. For 
example, when Dig = 591 = 1, Di4 = 0, 810 = 801 = 2 and 

= = x = 1 
811 l, we get p q dica 


Sii 


III. THE PROGRAM OF ROESSER'S EQUATIONS 
WITH SCALAR COEFFICIENTS (FIRST ORDER) 


A. AN EXAMPLE 
For a 4x4 data field the S and R matrices are indexed 


as follows: 





ο Πα Εμ R matrix 
For 4 x4 Matrices 
The Initial Conditions are given by the values 
R(1,1)4 ιδ ٧) 1 ال ال‎ 


S(l,l), Stl72) XC 


The 2-D state variable equations can be written as: 


R(i+l,j) = AjR(i,j) T A4S(i,j) سل‎ B,u(i,j) 
S(i,jtl) = AjR(i,]j) + A,S(1,j) + Bou(i,j) 
Lr 
γίνη) -. ος C, ] 
5117 


The input 2-D field is taken to be, 


u(i,j) 1, for i-2j^-1 


= 0, otherwise: 


32 


The output data field is indexed as: 





Y output matrix 


CC THE 2D FOURIER TRANSFORM 
ED diScrEEete EOUrIeE transform Y(m,n) of the output 


y(1,J) can be written as, 


M-1 N-1 -ἠ2π--τ ز‎ 
Y(m,n) = اا(‎ 6 ۰ 
2=0 k=0 
or for convenience, 
Μ N MUN Eo 
MGE) = ) ) y(£,k)e e 
Q=1 k=1 
ο 2-) D.F.T. {y(i,j)} 


MxN: The dimension of the given data vy(£,k) and D.F.T. 
ὑπ ο 21 


y(£,k): Given data (The output as described above). 
To develop the D.F.T. for two-dimensional signals we 


consider a finite area sequence y(%,k) which is zero outside 


دد 


the interval 0 < 2 < M-1, ο ο ο πο a Lu 


(M,N) and construct the periodic sequence: 
y(£,k) = YiC(2)) CK) ) I 


The original sequence y(%,k) 1S recovered by extracting one 


persod of ylk) دل‎ 


~r 


VIL y(2,k)R, e) 


QN 


ال n‏ 0 مث تت 3 ول 


RU Di Ξ 
Pun 0, otherwise 


We then define the discrete Fourier transform of y(Z£,k) to 
correspond to the Fourier series coefficients of ΤΙ 
However, just as we did with one-dimensional sequences, we 
will maintain the duality between the time and frequency 
domains by interpreting the D.F.T. coefficients to also be a 


finite 2-D sequence. Thus with Y(m,n) denoting the D.F.T. 


of y(2%,k), we can write 


٤-1 N-1 - 520 -j20 
πα... ) ) y(%,k)e e Ru ν Un, n) 
£20 k=0 í 
Or 
; M-l N-l jan 12111; 
VOLT = MN ) Ü Um) e ^.^ Ry xU IKE) 
m=0 y=0 
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Or, 


As an example, consider the case for M = N = 5. 


Given 2-D Data Sequence 


= 


2 ۹10 سه‎ 2 2 ASIA 2,5 Matrix 5 ۶ 
pu eK) = SINE. ٢ د 4 د د دن‎ M=5 N=5 
K IA, e AIA, 42 4:5 1-1 2 4, 51 د1‎ DS 


nen, 
οι πα (D O να - 2907 5) 
TA ceil: Views) ο τι ο Έτσι ay (Z, >) 
BUE, - I SENE τα 2( ۳ درو)‎ y 4) a y (3,5) 
πε] ,ΏΞ] 


+y(4,1) + y(4,2) + y(4,3) + y(4,4) + y(4,5) 


ΡΕ νο 2 στις ο ο ο) τν») 


B5 


L-)8 2(-)9 


(-)8 (js 


S 


—L-)8 (Π--) 9 


5 


---)8 (Π--) 9 


3(G‘T)A + 


κος e 


m 9 


Ə(p'S)Å 
o(pp)A 
ar EA 
۸م )ه‎ 


AE 


EU 


)1 (٨٩ 


(Tî) A+ 


(TE) A+ 


(Tz) A+ 


(I^D4 


(14S) A+ 


(T’7) At 


(I'£) A- 


(TZ) A+ 


(T’T)A 


£-u' au 


= (ETA 


نند "لام 


- ٨75 
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ς 2(G’SG6A + τ a(p’s)A + τ a(e’s)A + ς a(7‘S)A + s 9(I'S) A« 
(“C-) 8 (Π--)}8 (Π--}8 —E-)8 (7E-)8 
ς 9(S'p)A 1 ς o(p'p) A + ς o(£'v)A 4 s o(c'p)A )د‎ ς Ə(T“F)Z+ 
n 9 (---)9 (S (Π--}9 πο. 
ς 9θ(ς/ε)ά - ς 9(p'c)A + ς a(e’e)A + F Ə(Z'£)Ä + : 3(T’¢) At 
mr zL-)v CEST ye)? zL-)v 
ς Ə(G'Z)Å + ς 9(p'z)A * ς o(c'z)A * ς ϱ(σ’ζ) ἆ + τ (1'7) A+ 
m ^ πε-)ζ (qe eae E 
(GDA + (P'T)ZA + (ε/τ) ἆ + (Z’T)A + ) A 
ς Θ(ς/ς)ά - ς Ə(p'S)Å + ς θ(ς/ς)ἆ - τ Ə(zZ*G)Á + (T'S)Z+ 
={-) pz =C-) gT (=C-) ZT (=C-)9 
c Ə(S'p)Ä + τ Ə(p'p)Ä + ς a(e‘pA + د‎ Ə(Z'p)Ä + (t'v) 
-rz (=C-) 8T (Cozi (=C-)9 
ς a(g‘e)A + P a(p’e)A + : SE ENA τ a(z‘e)A + (TE) A+ 
(—C-) $2 (=C-) 8T (=C-)2T ={-)9 
: 9(g'z)À 1 š Θ(ρ}σ)ά 4 è ο(ς’ ζ)ἆ -- : (ج)) ې ۸(ح 'ج)ه‎ ٤ 
—-)vZ (=C-) gT (=C-)2T ={-)9 
ς Ə(S'T)Å + : 9(P'DA + τ a(e’T)A + τ απ“ 


=C-) #2 (-C-) 81 -(-)zI ={-)9 


οσα 


= (Τζ) 


p-u' Tau 


Φτ)λ‏ د 


>y 


In Appendix A we give a listing of the programs that have 


been written to generate y(i,j) and Y(m,n). 


C. NUMERICAL EXAMPLES 
Three numerical examples which depend on Equation (II.19) 
are used to demonstrate the program in Appendix A. 


First example: 


ο E.) 
H(Z,,25) = ποτ πα ου 
224 «32. 
yields 
a 79 
AJ = 291 = 0.3 
A, = a1 0 = 0.2 
Ci Ξ bio = 0.5 
C> = bol Ξ ο 
By = 1 
B. = 1 
D ERU 


After substitution of these values in Eqs. (II.16) and τπτ, 


we identify 
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A, ο = ὃς 

= q = 0.3 

uo - V‏ وا لا 

ou = bol = 065 ΜΙ ΠΠ 
q 519. - 5» 

By = 1 

Bo = 1 

D = 0 


Second example: 


Proceeding in a similar way with 


ccc cod MT 
H(z4;25) = ٥٢ ee - يپ حر‎ 
1-0.1252, 0.22. "0.171 z 
yields 
811 = 0.1 boo = 0 
AJ = 801 = 0.2 Dii = 0.2 
A, = p = 0.125 — Ὁ 3 
A, = q = I ος = 0.15 ΤΠ ΠΟ, 
A, = 810 τα. 1 5 
C4 = bol = 03 
C, = Dio = 0.25 


2د 


Third example: 


Proceeding in a similar way with 


01252 ο ο ٢ ο 
H(Z 22) ' =] : -I a = -1 
1-0.1352i -0.25z, -0.1524 z 
yields 
boo = 0 Dii = 0.25 
891 7 05 
A1 EE a O25 
A —-— pos (0R 
Bo v E 1.4 (ITT μα 
A, = ajo = O J 
οι 7 bol = 0.15 
5. D1p = 0.25 
B, = 1 
Bo = 1 
D = 0 


Zero initial conditions were assumed for all examples. The 


simulation results are presented in Figures 3-l, 3-2 and 3-3. 
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VATION: 
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URIGIN 


به 


for Example l 


Y (m,n) 


Sequence, 


Ομ 


Pigure la: 


t^ — 
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N x. ` τ N N Na. ` — ^ N 
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b N ` ` x NS δν "` — -^ | 
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Contour Map for Figure 3-la 


πο oe JID. 
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for Example 2 


Sequence, Y (m,n) 


ΙΙ. 


Figure 3-28. 
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MAP 


CON FEOUR 


Contour Map for Figure 3-2a 


Figure 3-2b. 
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Sequence, Y(m,n) for Example 3 
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Figure 3-3a. 
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lnEorderaPo verity the correctness of the output produced 
by Roesser, the 2-10 D.F.T. Y(m,n) plots for these examples 
were compared with the corresponding ΓΗ (21922) |. 2-D transfer 
Function plots ΓΗ (2122) | for Examples 1, 2 and 3 are shown 
in Figs. 3-4a,b, 3-5a,b and 3-6a,b respectively. The listing 
of a program used to generate these plots can be found in 


Appendix B. 
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for Example 1 
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Figure 3-4a. 
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3-4b. Contour Map for Figure 3-4a 
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Figure 3-5a. 
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Figure 3-6b. 


IV. EXTENSION OF ROESSER'S MODEL TO SECOND AND HIGHER ORDERS 


A. MINIMIZING THE NUMBER OF SHIFT OPERATORS 
In order to minimize the number of shift operators we follow 
the procedure given in Kung [Ref. 8]. Let us consider the 


Simple 2-D IIR filter transfer function given by 


e دی‎ 
11 1 52 2127717 
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TO L 012 2 JO I 2121 2 


zu =) 
b + Z. + Z 
00 P1021 Fo 
H(z , ) 
πο 
l-a 


2) l 
— لکت‎ ( LV ٢ 


B (Z, 2 


Our problem will be drawing a detailed signal flowgraph for 


the system function H(Z,,2.). We can do this simply enough by 


combining the flowgraphs on Figures 2-2 and 2-3 to get the 





Figure 4-1 


This flowgraph can be made even simpler because the shift opera- 
tion is distributive over addition. We can combine the two 22 
operators into a single one, yielding the following flowgraph. 
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OUTPUT 





Figure 4-2 
Doing so reduces the number of shift operators that need to be 
implemented and consequently the amount of storage necessary. 
There are other signal flowgraphs which E rise to the 
desired system function H(Z1;25). For example, we could 
wert the order of the B(z Z3) η σα the feedback 


kop Containing له اب د‎ to obtain the block diagram: 





Figure 4-3 


Then, when we substitute Figures 2-2 and 2-3 for the blocks as 
before, the two E chains will contain the same data and can 
be merged to yield the signal flowgraph in Figure 4-4. This 
glowgraph has a total of four shift operators, and it mini- 


mizes the number of an operators. 
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Figure 4-4 


zii 
1 
operators may be obtained from Figure 4-4by the 2-D transposition 


Another signal flowgraph that minimizes the number of z 


theorem to obtain a transposed network. Like tts 1-D counter- 
part [Ref. 9], the 2-D transposition theorem states that the 
transposed network, which is obtained by reversing the 
directions of all) Ehe arrows ima signal 1٧ ٢ لل‎ πα. 
the same system function as the original network. If we 
reverse the direction of all the arrows in Figure 4-4 and then 
redraw the flow graph with the input port on the left and the 


output port on the right, we get the flowgraph shown in Figure 4-5. 





Figure 4-5 
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This transposed flowgraph may be preferred in implementations 
with limited wordlengths since the attenuation due to the 
) 221085" که‎ H (21525) occurs before the gain due to the 
"poles" thus lessening somewhat the possibility of arithmetic 
overflow in the intermediate computations. 

Using the notion of transposition at both the flowgraph 
level and the block diagram level (note that Figure 2-2 is the 
transpose of Figure 2-1) the flowgraph can be manipulated to 


yield a realization that minimizes the total number of shift 


operators. 

As we Saw earlier, however, a zi operator will require 
substantially more storage than a Zi Operator for a sow py 
row ordering of input samples. Consequently, it may be more 


economical to minimize not the total number of shift operators 


2 


If the filter is realized by using a separate micro- 


(as in the 1-D case) but the number of z operators. 
processor to compute samples of.each node signal, storage 
may be less of an issue. 

In this case, we may want to minimize the total number of 
nodes in a flowgraph in order to reduce the number of micro- 
processors in an implementation. 

As digital technology progresses, the relative costs of 
storage, computation, and interconnectivity keep changing. In 
the future digital systems designers may have radically differ- 


emtecriteria for optimizing a filter realization. 
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B. A SECOND ORDER ΜΠΟΡΕΙΣ 
Looking at the flowgraph in Figure 4-5 and developing a 


state variable implementation from it, we shall call the output 


of the top 2 Operator R4 (1,3), the output of the lower 


- operator R, (i,j), the output of the left ES Operator 


S (i,j) and the Oütput of ene right zo operator S5 (1,3) as 


indicated: 
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In general the foregoing equations can be written as: 
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Now we can give an expanded version of (IV-2): 
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Equations (IV.2) and (IV.3) represent an algorithm for comput- 
ing the samples of the output signal from the samples of the 
input signal. Just as in the preceding subsection, the amount 
of memory required to store the state variables depends on the 
order in which the output samples are to be computed. It is 
possible to envision a m sss architecture for comput- 
ing equation (IV.4) by assigning each processor the responsi- 
bility of computing the next value of a particular state 
variable given the current input value and the current state- 
variable values. Equation (IV.3) could be implemented by a 
filter microprocessor to generate the desired output signal 
values. 

In such an architecture, m ana of the number of 
microprocessors corresponds to the minimization of the number 
of state variables, a problem studied thoroughly in the 
literature. Other state-variable forms with the same number 
of state variables can also be found that will realize the 
same system function H (21525) and may exhibit lower coefficients 
Of sensitivity or round-off noise [Refs. 2,10]. 

For the RES 3956 ۴2 210 205 2-D IR filters tnde 15, 


iil ers with a system function of the form: 


Όρο 


στ το. - ee 
l1 2 A(Z, +2.) 


where boo is a constant and A(Z, 125) is a 2-D polynomial, it 


can be shown that state variable realizations based on signal 
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flowgraphs, using the output of the shift operators as the 
state variables, require the minimum number of state variables. 
They are minimal realizations [Ref. 2]. 

From the above equations corresponding to the second order 
Roesser model, the program in Appendix C, was written. This 
program uses the values of coefficients of Η(Ζ1 22) as inputs 
and it generates an output, y(i,j). Next, the program finds 
the 2-D Fourier transform of this output matrix, and compares 


JE 


it to the transfer function H (04 ^04 
Numerical Example 4 
In the following three examples (first and second orders), 


we use the coefficients of first and second order transfer 


functions. We consider the special case of "all-pole" 2-D IIR 
filters, i.e., filters with a transfer function of the form: 
le 
00 l 
Ha aa ر‎ ae ست تت‎ 
1 : Ἀ(Ζ1/22) Ἀ(Ζ122) 

where boo is constant (unity in our case) and A(z, 12.) is a 
2-D polynomial. It can be shown that state variable realiza- 


tions based on signal flowgraphs, using the output of the shift 
operators as the state variables, require the minimum number of 
state variables. They are minimal realizations [Ref. 2]. 

For the third program we have a graph for the case of a 
BP filter. 


Example #4 


1 - 
pz = —S[V T vs Oo I— (1V:6) 
i 2 EI -I -I.-1 
1-0.224 -0.52 “0.171 Ζ2 


Η(Ζ 
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The |Y(m,n)| for this example is plotted in Fig. 4-la. The 
corresponding contour map of 2-D surface is shown in Figure 
Z5. 


Example #5 


= 1 
Η(Ζη 25) πο απ a a A a I CIV 


- τ. a T= 
1 0.0.7] 0.3452, 0.1252, Z 0.1z; Z 5 


miew2—-) D.F.T.\|y(m,n)| of the output of this filter is shown 
miele 4—2a. The corresponding Contour map is shown in 
Fig. 4-2b. Ν 


Example #6 


-0.12540.252]1«0.125221-0.1252] ^27. 
52) - τ 
a: ae 


= ae 
2 *0.12524 z 


NTS) 


pecan) | for this examle and the corresponding contour map are 
Shown in Figs. 4-3a and 4-3b, respectively. 

For reasons of verification, as before, Y(m,n) was com- 
pared to the actual transfer function H(w rw”) fOr examples 
4,5,6. These transfer functions plots and the corresponding 
contour maps are shown in Fig. 4-4a,b, Fig. 4-5a,b and Fig. 
4-6a,b for the examples 4,5,6, respectively. 

Gr TENSION OF THE 2-D STATE SPACE MODELS TO HIGHER ORDER 

TRANSFER FUNCTIONS 

IEEE Udtroduction 

During recent years, several authors (Attasi [Ref. ll, 


Fozmasimi and Mazchesini [Ref. 13], Givone and Roesser [Ref. 
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2-D D.F.T. Sequences Y(m,n) for Example 5 


Figure 4-2a. 
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Squences, Y(m,n) 
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Figure 4-3a. 
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Figure 4-3b. 
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17]) have proposed different state space models for 2-D 
systems. They have also suggested some extensions of the 
usual l-D notions of controllability; 09367۷37111 sane 
minimality to the 2-D case. 

However, these results are not quite satisfactory. 

They either lack motivation τους the stave Sp ΙΘΞΝ ο 
duced or the notion of state-space is improperly defined. In 
Chapter II we started with a comparison of all the current 
models based on a practical (circuit-oriented) point of 

view and on a proper definition of state. It is shown that 

the model of Roesser is the most satisfactory, in that it is 
also the most general Since the Attasi and Fozmasimi@eMazcho اا‎ 
models can be imbedded in the Givone and Roesser model. 

In Chapter ii we porn 11 πως ἡ 03( 07 0111 ال‎ 
between 1-D and 2-D systems is that in the 2-D case a global 
State (which preserves all past information) and a local state 
(which gives us the size of the recursions of the 2-D filter) 
can be introduced. 

2. Extension for 2-D Systems 

In [Ref. 14], Fozmasimi and Mazchesimi use the algebraic 
point of view of "Nerode" equivalence. In this framework, the 
state space arises from the factorization of the 2-D input/ 
output map.  Fozmasimi and Mazchesini were the first to realize 
that a major difference between l-D and 2-D systems is that we 
can introduce a global state and a local state in the 2-D case. 

The global state (which is of infinite dimensions, in 


general) preserves all the past information while the local 
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state gives us the size of the recursions to be performed at 
each step by the 2-D filter. However, Fozmasimi and Mazchesini 
271 2ا‎ κο explore fully the structure of the global state and 
its relation to the local state, so that the state space model 
they introduced is unsatisfactory in the sense that what they 
introduce as the state is really only a "partial state" (as 
defined by Wololich [Ref. 15] for 1-D systems). Indeed, this 
partial state does not obey a first-order difference equation 
(the notion of first order difference equation for linear sys- 
ΜΓ. partially ordered sets has been defined by Mullans and 
Elliot in [Ref. 16]) Attasi's model suffers from the same 
drawback as the Fozmasimi and Mazchesini one. 

On the other hand, Givone and Roesser [Refs. 17,18,1] 
have used a "circuit approach" to the problem of state space 
realization for 2-D systems. They present a model in which 
the local state is divided into a horizontal and a vertical 
state which are propagated, respectively, horizontally and 
vertically by first-order difference equations. From this 
point of view the global state appears as the boundary condi- 
tion necessary to propagate the state-space equations. 

However, Roesser did not provide much motivation for 
the introduction of such a model and seemed unaware of the 
iIMMMECIrcurt interpretation of their model since they were not 
able to implement an arbitrary 2-D transfer function, say 

(2122) 


στη 


i 


1 


Mitra et al gave an answer in [Ref. 19] by presenting 
an implementation method for 2-D transfer functions using delay 
elements 217 and 2pm We shall see below that this approach is 


consistent with Roesser's model. It is shown in [Ref. 8] that 


Roesser's model appears naturally as a way to describe the local 


state properties. For a (n,m) 2-D transfer function, 
n m ER 
b(z,,Z.) D pi ee 
Hz 2. = i 2 = UN LUN (IV.9) 
π (2122) n m ας : 
d NN 45471 22 
1Ξ0 10 


exhibits some canonical state-space forms (controllability, 


observability), which can also be written as, 


= ΜΕ 
) b. (z+ τ. 
_ ἆξο 
Η(ζ1/22) = n (IV: 10% 
e 
1=0 
Without loss of generality, we can assume agg 7 l and we 
denote 
Ξ τς τ 


Thus, using 1-D realization technique, Η(ζ1/22) of Eq. (τν ο 


can be used as shown below in Fig. 4-7. 
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Figure 4-7 


The realization is almost achieved: in addition to the 
n-horizontal delay elements, we need only m vertical delay 
elements to implement the feedback gains E I O TT 
and m other vertical delay elements to implement the readout 

gains (b, (z^), 0, | weenie Tims) the complete realization 
Shown in Fig. 4-8 requires only n+2m dynamic elements. This 
realization is a standard (canonical) one; its structure is very 
Simple and it involves only real gains. Note also that we 

need fewer dynamic elements than was suggested by the imple- 
mentations of [Ref. 19]. 


As mentioned in Section (b), circuit implementations 


: e πι 
with delay elements z and z are in a one-to-one corres- 


1 2 
pondence with state-space models of Roesser's type. The 
Suepues of the τ᾽ delays are the horizontal states and the 
uu 


eumsputs ot the z delays are the vertical states. Thus the 
implementation of the following figure can be transformed 


readily into the following state-space model. 
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D. PROGRAM AND EXAMPLES FOR ROESSER'S EQUATIONS USING 
KUNG'S MODEL 


This program (Appendix D) takes as initial conditions one 
horizontal state and two vertical states. The order of hori- 
zontal states is given by N and the order of the vertical 
states by M. 

We give two examples, one for N = 2 and M = 2 (Example 7) 
(two orders for horizontal states and 2 orders for vertical 
states) and one for N = 4 and M = 3 (Example 8) (four orders 
for horizontal states and three orders for vertical states). 


The first example is for a matrix 2 x2 and the second example 


۸4 1 
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E. NUMERICAL EXAMPLES FOR KUNG'S MODEL 
The following presents three examples. The first one 
corresponds to an "all-pole" 2-D low-pass filter. The second 
one is an "all-zero" 2-D band-pass filter (sin 4 sin Wo). 
I s bird One jis also a band-pass filter. All these examples 
are second order. The outputs of these examples are produced 
using Kung's [Ref. 8] state-space model. In this formulation, 
for a second order system, we require two horizontal states-- 
Pawo) and R20i,3) and four vertical states, S1(1) (i,j), 
πι, 521117} απὰ 52{2) (1,3). The program listing 
for implementing this model is given in Appendix D. 
Example #9 
The system parameters and the initial conditions chosen 
for this example are as listed in Table 4.1. The 2-D D.F.T. 
ly(m,n)| of the output sequence y(i,j) produced by the program 
ΠΠ mendix D is shown in Fig. 4-9a. The corresponding contour 
Map is Shown in Fig. 4-9b. 
Example #10 
The parameter coefficients and the initial conditions for 
this example are listed in Table 4.2. The 2-D D.F.T. 
sequence |Y(m,n)| for this example is illustrated in Fig. 4-10a, 
and Fig4.10b shows the associated contour map. 
Example £11 
The parameter coefficients and the initial conditions for 
this example are listed in Table 4.3. The 2-D D.F.T. sequence 
IY(m,n)| for this example are illustrated in Fig. 4-lla and 


Figure 4-llb shows the associated contour map. 
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Once again, to verify the correctness of our program, 
the D.F.T. |Y(m,n)| was compared to JH(w, 05) |. H (04 (0) 
and the corresponding contour maps are shown in Fig. 4-12a,b, 
Fig. 4-13a,b and Fig. 4-14a,b for examples 9, 10 and Il, 


respectively. 


F. SUMMARY OF PROGRAMS DEVELOPED 
The programs which have been written, cover the following 


orders based upon the different models. 


Appendix Order Model # of States 
A lst Roesser IL Hori bar NE J vertical 
(e And Roesser 2 horizontal, W vertical 
D Multi-order Kung ٨ horizontali 2v vertical 


In order to check the program listing, the same first 
order example was used on all programs.  Identical results 
were obtained. Similarly, identical second order examples 


were used in Programs C and D and produced identical outputs. 
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Figure 4-12b. 
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VS USE 9 Ρα PACKAGE 


De SS PACK 

"SSPACK" is a "state space system package," [Ref. 21] that 
is an interactive, state-of-the-art, software package for the 
analysis, design, and display of one-dimensional state-space 
systems. The work which follows adapts this program so that 
it can be used to produce 2-D data fields from state space 
formulations. A brief description of SSPACK follows. 

SSPACK is useful for a variety of applications in signal 
processing and control [Ref. 22]. The package consists of a 
supervisor which controls the operation of the software anda 
set of independent programs which communicate using disk files. 
The core of the package are the pre- and post-processors. The 
state-space pre-processor (SSPREP) program aids in preparing 
μυ τς τον the individual algorithm programs.  [Refs. 23,24] The 
state space post-processor (SSPOST) program displays and analyzes 
ΤΕ ἐεσπ the algorithms.  SSPREP prompts with a series 
Of questions in a menu format. 

SSPOST is an interactive command-drive processor. It is 
designed to help interpret the output of the various SSPACK 
algorithms, and display time histories: 

A is the Nx by Nx state transition matrix; 
HENCE ShnSNNx by Nu input transrtron matrix; 


C is the Nz by Nx measurement matrix; 
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D is the NZ by Nu feedtnru πο σας, 
W is the Nx by Nw process noise matrix; 
V is the Nz by Nv measurement noise matrix. 
The SSPACK works in multi-order form, using the transfer 


function of the l-D digital filter: 











ما٨۷‏ 6 اا 0-ا| 
FILTER‏ 


|-Ὁ (πν-εΝ سا‎ Ce 
φως ΟΝΤΑ Ou < FUT DATA 






Figure 5-1 


The present objective is to use SSPACK with a 2-D input data 
field and through the same transfer function, l-D digital 
filter, to accomplish 2-D output data field. 
B. DESIGN OF 2-D DIGITAL FILTERS USING I-D DIGITAL PITER 
STRUCTURES 
The idea of using two types of dynamic elements is not very 
abstract; it is very natural in delay-differential systems. 
However, before considering ltsepreactical avpolicatlons seo 
image systems, two remarks have to be made. The first is be- 
cause the "spatial" dynamic elements seem unimplementable, and 
we need to replace them by time-delay elements. Secondly, in 
order to have a finite order, we shall only RE a bounded 
frame system, i.e., we assume that the picture frame of interest 
is an MxN frame (with vertical width M and horizontal length 


N). Note that in order to use time delay elements, we need 


Ti) 


Pest FO Lind a wapo Code a 2-D spatial system into a 1-D 
(discrete time) system and vice versa. 
Thus we propose the following system, composed of three 
subsystems in series: 
٥۸ تت‎ τοτε πα πε Scan Generator codes the 2-D spatial input 


mee l-D time data according to the mapping function 





Ex) ΙΝ... ΞΘ Νοκ ας Ναι ce 
oS ML 
where M and N are relatively prime integers. For example, we 
Somerdqer a 2-D input data u(i,j): 
(0,0) Mo (0,2) τυ ial) 
(1,0) (1,1) Em 
wee) = 
ad QI UNI ,1) (N-1,M-1) 
Scanning 


π ο λε οις uit, 7) 1s scanned to produce u(*) as follows: 


meme — 9(0,0),(0,1) ,(0,2) ,.22,(0,M-1),(1,0),({1,1),..., 
eM) > (N-1,0)... Gia teM-)) 


eee —00 Ll, 2eeM, Mtl, Me2, εν. μι απ), t = 1M + JN 


For example, 


ΠΊΕ 





u(1,J) = 0 0 0 
0 0 0 
yields 
y(t) = [1 ο Ὁ Ὁ ο νο ο 
ii) A l-D (discrete time) digital filter processes the 


l-D data generated. This subsystem is implemented by replacing 


zi by D z, by A in a 2cD ëireunit realiza rom εσ. UD 
controller form). ô and A are chosen as: 
M 
6 = Ὁ = M-units delay element 
A = p" = N-units delay element 


iii) The Output Frame Generator decodes the 1-D (discrete 
time) output of the 1-D digital filter described above into 
a 2-D (discrete-spatial) picture according to the inverse 


mapping of (V.1). 


(1(t),3(t)) = PE Mod IE P E T) )۷ 


where P is a unique integer such that PM-PN = l and 0 < P < N. 


This formula is given in [Ref. 2]. Alternately, we can compute 
( 9) as 

1 = tea 7 
and 

J = OUotient (E/N) 
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- 


For example we suppose t = 19 with N = 10 and M = 9. 


Iucccorresponuding value in the 2-D case will be i - Remainder {2Z} - 
= 9 and j = Quotient (το = 1. So in the 2-D case we will 
ο ο (i,j) = (9,1). 


Another Example: For M = 4 and N = 5, the single index £ 


will be mapped into (i,j) as: 





Macmorocedure fon implementing 2-D filters using l-D filter 


EEpuectures 1s as shown below in Fig. 5-2. 
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1-0 (Time) vS GRE) 
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2-ὺ iNeUT DATA 
F\ ELD 


Figure 3-2 


Γι» 


The index scanning is required for the input data so SSPACK 
can be carried out simply because the input is assumed to be 
a 2-D unit pulse. This is followed by implementing the 
corresponding 1-D filter of Fig. 5-3 using SSPACK to convert 
the l-D output from SSPACK to a program for output index 
mapping--written as shown in Appendix €. The 2-D Fourier 
transform of the resulting 2-D field is then computed. 
Considering a bounded frame (M xN} system it is interesting 
to know the dimension of the global state (or initial condi- 
tions) needed to process the M xN future data field. Since 
vertical states convey information vertically, all the verti- 
cal states along the X-axis are necessary initial conditions 
and their dimension is mN. Similarly, all the horizontal states 
along the Y-axis are necessary initial conditions (with 
dimension nM). They convey information horizontally. 
Therefore, in the bounded frame case a total number of 
mN+nM are needed to summarize the "past" information. This 
very same idea can be used again from a computational point 
of view. Indeed, the number of required storage elements for 
recursive computations is also equal to mN+mN if initial condi- 
tions are not zero. However, it is quite often the case that 
the system starts with zero initial conditions; the size of 
storage required is reduced to mN (respectively, nM) which is 
used to store the updated data row by row (respectively, column 
by column). No storage is needed for the rest of the initial 
conditions--nM horizontal states (respectively, mN vertical 


states) since they are assumed to be zero. This is consistent 


114 


ε-ς SINO 








Lay? 3 D Dh ` Ux) "S CX)'s S (2) NS 


T i ΠΚ 1 \ 
beer : COH FREI A ΓΩ 


oA) 
ways و‎ L2 ἔσἔ.;-; 
say "S 


Xy. 


PES 


with the results of Read [Ref. 24] derived from a direct 
polynomial approach. 

Another interesting observation concerns the dimension 
Of MED figital filter contained by our 2=D digital filter 
design discussed above. Since it needs nM unit-delays and 
mN-unit delays, the corresponding 1-D state-space also has a 
dimension equal to nM+mN. Note that, despite the high dimen- 
Sion of the corresponding 1-D filter, its high sparsity is very 
encouraging for further studies. In short, following the 
above method of designing a 2-D filter, for the first order 


case, 


4 1 
H(zi:25) د چم بوس شه‎ ccc a au 


sas z s 2] 2; 

Using the above approach we get the 1-D filter realization for 
this 2-D filter which turns out to be as shown in Fig. 5-3. 

The detailed matrix equations for realizing Eq. (V.3) 
using SSPACK can be written as, The SSPACK produces a l-D 
sequence, which converted into a 2-D sequence using the output 
index mapping formulae discussed earlier. The listing of a 
program which does this mapping is shown in Appendix E. 

After obtaining the valid 2-D output data sequence yíi,]j) 
we next compute its 2-D D.F.T. to produce |[Y(m,n)| which for this 
example is plotted in Fig. 5-4a. The corresponding contour 


map is as shown in Fig. 5-4b. 
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Contour Map for Figure 5-4a 


Figure 5-4b. 


Vp OCF SIONS 


This thesis has dealt with the problem of modelling 2-D 
data fields in the state-space domain. First of all we have 
pointed out the main problems associated with the extension of 
l-D time-discrete state-space models to 2-D data fields. The 
remaining part of the thesis has been divided primarily in 3 
Parts. 

In the first part we describe Roesser's [Ref. 5] approach to 
modelling 2-D systems in the state space domain. Extensive com- 
puter simulation results are presented to verify.the functioning 
of this approach. This modelling approach has been tried out 
for the scalar (1x1) as well as for higher order (2 x2) etc., 

2-D systems. 

The second part deals with a modification of Roesser's approach 
as described by Kung [Ref. 7]. The main advantage of this approach 
is that the 2-D state-space model can be realized as a 2-D 
ο ασ , σι or tantly, ΕΠ15 2-5 circuit realization can be 
implemented as a 1-D digital filter. Computer simultation studies 
that have been carried out substantiate the making of this model. 
mesi D filter realization Obtained in this part turns out to 
be a very convenient starting point for the nezt part of our 
effort, dealing with the use of the 1-D SSPACK commercial soft- 
ware package designed for dynamic system simulation. 

In the final part of the thesis, we make use of the 1-D 


filter realization of 2-D state-space model obtained in the second 


ie 


part, and implement this filter using SSPACK. Some additional 
programming effort reuigred for input and output mapping was 
necessary. Programs for converting 2-D input and output se- 
quences to 1-D have been written separately. In this 

fashion we have succeeded in extending the applicability of 
the SSPACK to simulating 2-D linear systems as well. Once 
again, detailed computer simulations have been carried out 


to verify the functioning of this modification of the SSPACK. 
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APPENDIX A 


2 ونج د 
٣ -- 5‏ 
ον = Z‏ 
D —iróes 1 ry Microsofti FORTRAN7Z7 Vs. 20 2/64‏ 
SS TGRAGE: =‏ 1 
SLARGE‏ 2 
ان i‏ سف خد چب خي خي ÀJ» JE JE JE te a ΒΕ -ᾱ- Ε΄ ΞΕ ΞΕ ΞΕ -ᾱ- E‏ چد ځد JE‏ ځد a JE E ΠΕ ΠΕ ἈΞ ΠΕ ΕΕ ΓΕ‏ د i i HEP‏ د ځد د ځد چد ΠΕ ΞΕ‏ ځد xe‏ ځد E JE dE Ab dE dE‏ خد خد ځه p ἄρ‏ 4 
* * = 
EN * THE PURPOSE CF THIS PROGRAM IS TO COMPUTE AND GRAPH THE ^‏ 
uot + EQUATIONS OF ROBERT P. ROESSER IN THE "DISCRETE STATE-SPECE *‏ 
Su * MODEL FOR LINEAR IMAGE PROCESSING". *‏ 
CIE * *‏ 
HORE * EVANGELOS THEOFILOU *‏ 
ځد ځچي خد ځد JE‏ ځي کل ځې ځف AP JP J‏ خد خد خي خد خد خد خد خد E te HE HE JE E JP HP‏ ځد د i JE b E JE dE ME JE P JP E E ΒΡΕ ΞΡ i‏ د HF FF HF‏ بد پد چ Je JE Jb JE Jb JP de Ae Hb 3» AP ie‏ نا 1 1 
E E PROGRAM 2D-DaTAa-riELD‏ 
چ to‏ 
VARIABLE DECLARATIONS «ο‏ +++ 8 14 
τι ο οἱ).‏ 2ل ورک) له REAL τι: ο. ο το)‏ = 
τος. ρα στο ο) VERIFEXCIS)!ZUZVCSIT)‏ * 16 
NTEGER MASKKEOOCONMEDTG CE DMEWNGT (E19‏ 17 
CHARACTER#1L ANSWER š‏ 18 
CHARACTER#20 CTEXT‏ 13 
20 
٧۱٩ م٢‎ δι =7, YURR/ 7.07,‏ ویرت mI DATA Κα ο ο τι ος ο‏ 
2e * LOW/1.9E3S/, IPROJ/9/,NRNG/100/‏ 
25 
C καθε η I N P RO GRAM FEE 3 4 3 38‏ 24 
ao‏ 
τ; 4-343 OSH THE RERUIRED VALUES FOR THE MODEL «.ε-ᾱ‏ 
τ... Εν ται VALUES ΕΠΗ THe FOLLOWING VARIABLES(s. #,..) 2°‏ کد Sm τῷ απ‏ 
EE TA‏ 9 1 له 29 
READ (+, =) Ai‏ = 
πρ ο. του “Aas 7‏ 2 
Sl READ (*.*) Ae‏ 
Te ARE OE OO I) n S |‏ 
y RERD (+, +) AZ‏ 
WRITE (FF, 99) A4: °‏ ڪت 
READ (#, #) AS‏ 25 
συν κε 233) οι *‏ 6« 
ST READ (+, +) Bl‏ 
ΤΕ OE SSI) O‏ 28 
m READ (#,*) Be‏ 
ARLE ool) “ει‏ 40 
READ (+, +) C1‏ 41 
σι ια κ 99) 1-5: '‏ 42 
ReaD ο we) Ce‏ ت 4 
S WRITE (#, #02)‏ 44 
e READ (+, +) N‏ 
TECON 512-259 GOTO S‏ 46 
47 
WRITE (#,211) ENTER ',N,’ INITIAL CONDITIONS FOR MATRIX R(#.#)'‏ 48 
DO 33 I = L. ñN‏ 43 
SO WRITE (#,403) °R(1,°,I,’) 3s ° :‏ 1 
al READ (*,*) R(1, 1)‏ 1 
CONTINUE‏ 93 = = 
ك 
z^ ARIE COE SIND "ENTERO OUUNSCOEINTIISESCONDITIONSSEOR MATRIX S(#. 4)’‏ 
DO 100 I = 1,N‏ = 
τν Ὁ‏ ولا URIE (E ον St‏ 3 1 
S RERD (κ το S(L, 1)‏ 1 
CONTINUE‏ 200 32 1 
ca‏ 
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61 
62 
6-5 
64 
65 
66 
e 
68 
693 
7ο 
71 
τς 
73 
74 
τα 
76 
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7 Micrasoft 2۶201۴۹۸۷77 (ات .تا‎ Ὁ 
WRITE (#*, 419) 
READ (#, 200) ANSWER 
D ((ΒΝΘΗΕΠ ο νυ) SER: {ANSWER "Eb. yD) cote. 


U = 1.9 
**++*+ COMPUTE R AND S MATRICES ***** 
DO 101 I 51. 
DO 191 J = 1,N 
IF (I*1 .LE. N) THEN 
R(I+1,J) = P1*R(I,J) + P2Ə*S(I1,J) + Bl*U 


ENDIF 
IF JS Ε.Ε. Ὁ THEN 
Stl, J+1) = ھت‎ O T ID EB ARS OD e Box 
ENDIF 
U = O. O 
CONTINUE 


«κκ FILL 0% s THE TWO DIMENTIONAL GRID OF CONTROL POINTS **=*+ 
DO 1DG Lo us l 
DG 102 3 دول ګت‎ 
Z(1,J) =O. O 
CONTINUE 


33333 COMPUTE Z MAOTRIX s # 
DO 103 I ت‎ 1 
Dg 10.5. ο. N 
2) ( Ci*R( SJ, I) + C2<S(J, 1) 
CONTINUE 


κε OUTPUT THE Z MATRIX «κκ 

WRITE (*, 205) "'x*sewo ὅ mMoB TOROLOX OMNEM Je 3e 3C 36 3€? 
WRITE (*, 213) 

DO. 10401 = aN 


WRITE (*,200) (Z(I,J), J 5 1,ND 
WRITE (*, 210) 
CONTINUE 


WRITE (#, 213) 


WRITE (#, 418) 
READ (*, 200) ANSWER 
IF <¢(CANSWER .NE. °Y") LAND. (ANSWER NE IID οἱ. -- 


αν ASK THE PARAMETERS FOR THE GRAPH ****+ 
WRITE (#7210) 

WRITE (+, +) +++» ENTER PLOT PARAMETERS ام‎ 
WRITE (#*, 405) 

READ (*, *) AZIM 

WRITE (#*, 406) 

READ (A, o) Εαν 

WRITE (#, 408) 

READ (*, *) ITRIm 

WRITE (*, 409) 

READ (3, *) IDIV 

WRITE (#, 411) 

READ Ce, 133) ΕΠΕ 

WRITE (#, 401) 

READ (30200) ANSWER 
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tT 
ΤΕ ΕΞ: Ὁ, ΠΠ. ANSWER MEST vy 
CALL PLOTS Qo 
EISE 
CALL PLOTS(O, 53, 33) 
ENDIF 


ΞΓΒΡΗΙΠΕΝΡΣΨΟ, ΕΠ 
n 


) THEN 


CALL WINDOW(XLOL, YLOL, XUAR, YUPR) 


*+*+++ DRAW THE MESH SURFACE OF THE GRAPH «+44 
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606 به‎ 1۳0٩ )6 2 6 ك .9 و2‎ 2-D DATA. FIELD*® ,O. ὃς 14) 


3334 OUTPUT THE GRAPH scs 

ΕΞ: 2). 9 ο ο 9399) 

WRITE (*,4142) 

RERD (4,200) ANSWER 

IF DGUGNSUEN C πο νυ) ΕΠ. ENSUES «τα ντο BOE 15 


WRITE (*, 417) . 
RERD(*,200) ANSWER 
IF ((ANSWER .EQ. °Y’) «08. (ANSWER .ED. 'y')) THEN 
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cee τοι ο 5 ΤΗΕ WO DIMENTIONASE GRID O* CONTROL POINTS sess 


DO LUG 51. = 1,31 
DO 08 J = ول‎ πι 
Zetia) Ὁ 
CONTINUE 
ΖΕΜΩΧ zs -3.^3Ezo 
ZEMIN = 3. QEZ 
DN = (N-1)/2.0 
pP = 6.285185 
ΤΕ 107 1 = 1,N 
DBE ον د‎ πμ 
RLEART = 0.0 
IMGPART DT 
DO 108 L = 1,N 
DO 108 K = 1,N 


MICID) = EOS SBR CSI) )عم‎ L=-DN=1)/N) 

Rl (2) 2] SINC—=e4 (E>. E CTS DNS DZD) 

πως لت‎ οσο οσον —-DN= 1) 7 N) 

RECS) —OETNG-B3sCOSLDGECJ-SDN-—IJ/N) 

RLPART = RLPART + Z(L,K)*(R1(1)*RC (1) 
=R CSI DRE CS) S 

IMGPART = IMGPART + Z(L, MK) Æ(Ri(1)#RŠZ(5) 
رح دد‎ κας) 

CONTINUE 


ZF(I,J) = SGRTC(RLPART «#2 + IMGPART #+#2) 
DE τετ ο ο οὐ ο ο etc 
ZEMAN = ZE 
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EN 
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ΕΝ 
CONTIN 


* 33€ 9€ 
WRITE (* 
WRITE (* 
DG 5 ل0‎ 

WRITE 
WRITE 
CONTINUE 
WRITE (* 


μπιτ κε 
READ (* 
IF (¢ (AN 


HHH 
WRITE 
WRITE 
WRITE 
READ 
WRITE 
READ 
WRITE 
READ 
WRITE 
READ 
WRITE 
READ 
WRITE 
READ 


dede t € ¢ 
18 πι 
CALL 
EISE 
CALL 
ENDIF 


WRITE 
READ 


CALL W 


DE ٢ 
ان‎ 
CALL 
DO 1 

LD 
LW 
CONT 
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Microsort5 FORTRANZT7 VS. Ze 


DIF 
(ΖΕ(1.3} LT. ZEMIN THEN 
ZFMIN = ZF (I, J) 
DIF 
UE 


OUTPUT THE ZF MATRIX +H 


COS) καπ FOURIER TRANSFORMATION ?",N," 
Scd 

= 1,N 

(*. 200) “(ZF Cle. J), ریم نت لل‎ 

(#210) 


813) 

, 418) 

,200) ANSWER 

SWER «ΝΕ. Y!) .AND. (ANSWER «ΝΕ. 7y?)) 


ASK THE PARAMETERS FOR THE GRAPH desse 


Jaga - 
“s 
:شتو‎ 
Q Oz/34 

۸ NM. o eee? 


00 16 


#E109) 

(η) Tee ENTER PLOT PARAMETERS φον) 
(#, 405) 

(w, *) AZIM 

(*, 406) 

(x, *) SLEY 

(*, 408) 

(I) ITRIM 

(*, 409) 

(+, +) IDIV 
(*,41i) ; 
(*. 129) MGE 
GE, 401) 

(*, 200) ANSWER 


INITIALIZE PLOT38 «κκ 

ANSWER EQ. °Y’) .OR. (ANSWER .EQ., °'y?)) 
PLOTS (9, 0, Z) 

BPLOTS(€O. 99. 999 

(*, 420) 

(*, 200) ANSWER 


INDOW (XLOL, YLOL, XUPR, YUPR) 


ANSWER ER. 'Y') .QR. (ANSWER .EQ. "'y?*)) 


= (ZFMAX—-ZFMIN) /FLOAT (N) 
ZLEVEL(ZP, 31, 31,N,N, OLeY, ZLEV, Mel? 
10 1 = 1,N+1 
Is TM 
GT(I) = 1 
INUE 


* DRAW THE CONTOUR MAD wwwx 


THEN 


TOEN 


CALL ZCNTUR(ZF,31,31,N.N,0.5, 0.5, 7.5, 5. 5, ZLEV, LDIG, LWGT, 


CALL 


Ney Onl One 
SYMEGL<S. 5S, 0.0, 0. 2. CONTOUR, dee د‎ 
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SIM 


Name 


Al 


00 
OOM £ (C) fu ^ Cu O NYO (CO FF CL 


[ú [O D3 fO fo fe D (ó Γι fo fo fo fo fu 
m m m OQ O O cn cn c cn ca cn 


(n 
“ 


3 
2693 
270 
= 71 
= 
273 
274 
m 
EG 
277 
278 
τση 
280 
231 
cc 
29S 
=84 
eB 
= 56 
EO 
288 
Ed 
290 


231 


16 


1 
ZOO 
zog 
210 
ت‎ 
els 
213 
300 
Seog 
400 
401 
402 
403 
404 
405 
406 
408 
409 
411 
412 
41 
417 
418 
4193 
429 


Type 


RE SL 


1 a 
2 ٩٢1 اځ نه ه٣ ح‎ FORTRGN77 VZ. X0 oz 


ELSE 
333-3 DRAW THE MESH SURFACE OF THE GRAPH »**«- 
ت8‎ 5 ١) و2 2 .9 و22 .اا وم 2 و1 2 ,اا ول ول وات لرک‎ SS. SS, DINO, 
S [PROJ, 1, ZLOW, =, ITRIM, MASK, VERTEX) 
κ ANNOTATION CF THE GRAPH xe # 
EAE SYMEOL Gara, Oss. Oe, AZIMUTH: ο ο τω) 
Εν Εαν του οσο ο ο Ε "σι ο, ος 5) 
Cmuier SV MMOL ao. a, O.0,0..2, ELEVATION: ο. ος 1) 
EAE ٧816-2226 2339/۱9 و2‎ 25۱١ 09 OED 
DY = (ZE(1, 1)/30.9) + ELEV 
)لت‎ DEADIN O I O ZE CI DEM, AR, YR) 
ΕΕΣ SYMBOL ο το ο. τος -'. 0, 0. L) 
EAE SM EOE I Oo, O71,0. 2, * = ORIGIN’ 70.0, 10) 
ENDIF 


CALL SYMBOL (1.0,6.75, 0.25 
CALL SYMBOL (6.0,6. 5,0. 2," 


Sea Oe ie EOD 


7 
ΠΕΡ": 


x» OUTPUT THE GRAPH ***++ 

SE MOG O ο 999) 

WRITE (*, 413) 

READ  (*,200) ANSWER 

IF ((ANSWER .EG. 'Y') .OR. (ANSWER .EG. *y')» GOTO 50 
ENDIF 
WRITE (*, 413) 
READ (*,200) ANSWER 
ο ο... »-6Q, 77?) ο. “ANSWER .EG. *°y’)) GOTO t9 
STOP 


FORMAT (AZO) 

FORMAT (A) 

ΕΟΗΜΑΩΤ{/, Ξ0Χ, A25, I2, aZ, 12,8, /) 

FORMAT t /) 

FORMAT (/, AB, 12, A47) 

FORMAT (/, 2X,’ (AZIMUTH 320.0)", 46X,’ (AZIMUTH 230.0)°, /) 
FORMAT (/,2X,* (AZIMUTH OSO. 0) ' , 46X, ° (AZIMUTH 140.0)°,7) 
FORMAT (10(F7. 2, 1X)) 

FORMAT (/, SX, A4, \) 

FORMAT (9X, \) 

FORMAT(/,SX,'SEND GRAPH TO THE PRINTER(Y or N): ',N) 
FORMAT (/, SX, ' NUMBER OF ROWS/COLUMNS FOR R AND S(1 to 25): °,\) 
FORMATS. Aá, I2, 03, N) 

FORMAT (SX, AZ, IZ, AS, \) 

FORMAT (/, 5X,’ AZIMUTH(O.@ to 360.0 DEGREES): ',N) 
FORMAT (/,SX, 'ELEVATION(30.0 to -30.0 DEGREES): ',N) 
FORMAT (/, SX,’ TRIM(O=NO, 1=X%s, Z£2Ys) : 7, \) 

FORMAT (/,5X,’2,4 OR 8 SUBGRIDS: °,\) 

FORMAT (/,5X,’° TITLE OF GRAPH(UP TO 20 CHAR): °,\) 
FORMAT(/,SX,'DO YOU WANT TO CHANGE PARAMETERS? 7°, \) 
FORMAT(/,SX,'DO YOU WANT TO REPEAT THE PROCESS? °, \) 
FORMAT(/,SEX,'DO YOU WANT FOURIER TRANSFORMATION ? °,\) 
FORMAT (/,5X,’DO YOU WANT TO MAKE GRAPH ? '*,N) 
FORMAT(/,SX,’°DO YOU WANT TO CORRECT ? ',N) 

FORMAT (/,5X,’DO YOU WANT CONTOUR MAP ? ,Ν) 

END 


Offset P Class 
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D Line# 1 7 
Az REAL 

As REAL 

A4 REAL 
ANSWER CHAR*1 
AZIM REAL 

ἘΠῚ REAL 

Be REAL 

T RERL 

Cz REAL 

ES 

CTEXT CHAR*#EO 
DLEV REAL 

DN REAL 

DY REAL 

251 REAL 
FLOAT 

I INTEGER*cz 
[Diy INTEGER*Z 
IMGPAR REAL 
IPROJ INTEGER*¥2 
ITRIM εσεις 
ال‎ INTEGER*2 
K INTEGER*2 
Ὲ INTEGER*¥2 
LDIG INTEGERS 
LWGT INTEGER? 
MASK INTEGER#2 
N INTEGER*2 | 
NRNG INTEGER*¥2 
P RERL 

R REAL 

R1 REAL 

R2 REAL 
RLPART REAL 

S RERL 

SIN 

SQRT 

REAL‏ لا 
VERTEX REAL‏ 

χι. REAL 

XR REAL 

XUPR REAL 

YEDE REAL 

YR REAL 

YUPR REAL 

Z RERL 

ZF REAL 
ZFMAX RERL 
ZEMIN REAL 

ZLEV REAL 

ZLOW REAL 

Name Type 
MAIN 

MESHS 

NUMBER 


50 
34 
so 
74 
114 
42 
46 
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با 


g^ 


126 
erg 
166 
146 
118 


60 
124 
130 

ee 
5 

86 
ede 
134 

6000 
6062 
0 

ΞΘ 

24 
179 

O 
0 
8 
186 


2500 


76 

Ó 

e 
1S0 
10 

6 
154 
14 
3000 
8844 
138 
ντ 
12688 
18 


DNEDTRINSTIS 


INTRINSIC 


LARGE 
LARGE 
LARGE 


LARGE 
LARGE 
LARGE 
LARGE 
INTRINSIC 
INTRINSIC 


LARGE 


LARGE 
LARGE 


LARGE 


Class 


PROGRAM 


SUBROUTINE 
SUBRGUTING 


126 
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SYMBCL 

WINDOW 

ZENTUR 
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Pass One 
2391 Source Lines 


A> 


SUBROUTINE 
SUBROUTINE 
SUBROUTINE 
SUBROUTINE 
SUBROUTINE 
SUBROUTINE 
SUBROUTINE 


No Errors Detected 
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rm 
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(η 
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(o 
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[u 
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eo 
24 
im e 
GE 
2m 
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a 
DE 
=4 
ied 
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Sr 
D8 
39 
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41 
ac 
4 2 
44 
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SO 
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STORAGE: = 
BCFAGESIZE: 53 
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THE PURPOSE OF THIS PROGRAM 15 TO COMPUTE AND GRAPH THE 
FREQUENCY RESPONSE J^" sa ο ο TEE. 


ok ok‏ ې ې 


EVANGELOS THEOFILOU 
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PROGRAM 2D-DATA-FIELD 


Heee YORIABLE DECLARATIONS +++ 

REAL AUT TIBI RIOT ثق‎ Sy Ro e n s 
RLPART, IMGPART, Z2(5i,51), 
ΕΠΕ. ο ED 

INTEGER MASK (ZOOO), LDIG(51) , Ewer (51) 

CHARACTER æ+1 ANSWER 

CHARACTER«zO CTEXT 

Dara χι οι ο  Ὁ το و927‎ 1227 8 πι αυ”. 

ZLOW/ 1. OESS/, IPROJ/9O/,NRNG/100/ 


η I OON DO ROG R Am ssa d 1F 3 1F 3 1 1F‏ په نپ 4F‏ یاځياځیا څر 


WRITE (#,401) 

READ (*, *) Ef 

ie CIT CG. ES) GOM ο 

WRITE (#, 402) 

READ (*,*) K 

K = K +* 1 

WRITE(*,*) ? 

DO 100 1 909041 
DO 100-.J πο m 


ENTER VALUES Or COEXIGIENTSI 


WRITE(C(CM ADAMS I. v. T T 
READ (*, *) B(I+1,J+1) 
CONTINUE 


لح وه سح لل 191 DO‏ 
e KS‏ اا DO‏ 


WRITE (A لور ولوت 2 0457د‎ 
READ (+, *#) ACI+#1, J+1) 
CONTINUE 


WRITE (*, 419) 
READ  (*,200) ANSWER 
IF ((QNSWER .EQ. *'Y'!) .OR. (ANSWER .€0. 'y*')) GOTO το 


A FILL O's THE TWO DIMENTIONAGL GRID OF CONTROL POINTS s**-* 
DG 15027 r= pis 
DO 1192 3 اولح‎ 
ζει ου ο. 
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JE 3 aP te 
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= 
- 


60 
El 
62 


- 
-— 


64 


= 
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66 
E 
ea 
63 
70 
real 
ne 


- 
-— 


74 


e 
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76 
r 
793 
2 
80 
81 
az 


< 
—_ 


34 


= 
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86 
87 
88 
5 2 
20 
21 


— 
-= 


as 
34 


= 
= 


36 
2m 
ES 
22 
100 
101 
KOS 


107 


104 


10:3 
102 


5-58 z 
Em 
21309326 
7 MIE SRE 2 27 7 7 νο 0-7 5 
CONTINUE 
ZMIN = 3. رت258‎ 
ZMAX  -8. ΞΕΞΟ 
2 = πο 772 
STEP = 2«D / εἶτ-ι) 
Wi = -0 - STEP 
L πο 
πο 12 e i. ID 
Wi = Wl + STEP 
We = -ὂ - STEP 
DO IOS TO SIT 
ټ با نت با‎ 8# 
We = We + STEP 
18 O41 = OL A=s 
DO 104 N = O,=Z-1 
RI (M+1, N+1, 1) = CIS(-M * Wl) 
ALAMEIN 2 SIN(-M + W1) i 
2(Me1,N+1,1) = COS(-N +* We) 
R2(M+1,N+1,2) = SIN(-N * WE) 
CONTINUE 
RLNOM = 0.0 
IMGNOM = 0.9 
ALDEN = 9.0 
IMGDEN = 0.0 
DO 105 M = 0,K-1 
DOMOS N = 0, K-i 
RLNOM = RLNOM+B(M4+1,N+i) #(R1LOM+41, Nei, 1) *#R2(M+1,N+i, 1 
— Ri(M-1,N-—1,2) *RZ(M-1,N-1,2)) 
IMGNOM = IMGNOM+E (mM+1, N41) * (RL (mM+1, N41, 2) #ROCM+1, N41, Z) 
αν οι Ότι RL Med NS 1) 
RLDEN = RLDEN+A (M41, Nel) *#(R1 CM+1, Nel, 1) #RECM+41, Nel, 1) 
چا‎ κα ) ۳٣و‎ ٩٣ ول )2 بل« زګ‎ ل١‎ =) 
IMGDEN = IMGDEN+A (M41, Nel) خم‎ (RL CMe, Nel, 2) *#RECM41, Nel, =) 
+ RƏ(M+1,N+1, 1) *R1 (N+1 , N+1, Š) ) 
CONTINUE 


ELEMENT - SQRT(RLNOY**c + IMGN 


SORT (RLDEN««Z 
Z(I,J) = ELEMENT 
τι: Πακ THEN 
ZMAX = Z(I, 
SNDIZ 
IC Z0 ο το. EN 
EMIN σοι 
ENDIF 
CONTINUE 
CONTINUE 
«κκ DUTDUT THE Z MATRIX sew 
MAITE (x, 20S) sss: 27 m Q T R +: 
WRITE (*,2192) 


Pol 


QM33*2) / 


+ IMGDEN»**Z) 
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7 Microsoft 7ORTRAN77 Va. 
DO 106 I - 1,IT 
WRITE Ct, SOU) σα ντ لل‎ o, 
WRITE (#*, 210) 
CONTINUE 
WRITS (*,212) 


WRITE (*, 418) 


READ  (*,200) ANSWER 
IF ((ANSWER .NE. 'Y’) J AND. (ANSWER .NE. 'y')) GOTO 15 


κκ ASK THE PARAMETERS FOR THE GRAPH λενε 


WRITE (*, 210) 

WRITE (+, æ+) |? +++ EN TER PLOT PARAMETERS #4? 
WRITE (*, 4109) 

READ (*, *) AZIM 

WRITE (*, 411) 

READ (*, (عد‎ EREY 
WRITE (*, 413) 

READ (*, *) ITRIM 

WRITE (*, 414) 

READ (*, κ) IDIV 

WRITE (*, 415) 

READ (*, 129) ETEXT 

WRITE (*, 51) 

READ fg AOD ANSWER 


weeee INITIALIZE PLOTS38 s 

167 Ὁαβδίεα ἙΈσ YN One ΟΕΕ, εν, ,πεὴ 
CALL PLOTS(O, 0, &) 

ELSE 
GALE PLOTS (to, 33, 52) 

ENDIF 


WRITE (*, 420) 
READ (*,200) ANSWER 


CALL WINDOW(XLOL, YLOL, XUPR, YUPR) 


(ANSWER .EQ. ?Y?) .OR. (ANSWER .EQ. "v )) THEN‏ ( م1 
οι Εν = (ZMAX-ZMIN) /FLOAT (IT)‏ 
GRLL ZLEVEL(Z, 91, S91, IT, IT.DELEV τι.‏ 
DO 108 I = 1, IT+1‏ 
EDIGSCI) =‏ 
LWGT (I) 1‏ 
CONTINUE‏ 
DRAW THE CONTOUR MAD s #e‏ 3333€ 
CALL ZCNTUR(XZ, 591, 91, IT,.IT;0.59, 04d 259 SS SR ASE EBTER ENOT‏ 
ο, το του‏ 11-1 
CONTOUR ΕΟ...‏ وت .لا ولا .اوه CALL SYMBOL CS.‏ 
zm‏ 


Jo 


(Tı al, (-1 


hy 
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153 
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ΙΕ. 
τος 
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164 
156. 
166 
Do 
168 
152 
170 
D 
Dc 
τς 
174 
175 
176 
TET 
178 
793 
180 
E21 
195 
18. 
194 
185 
186 
Το; 
188 
se = 
1:90) 
121 
LEE 
193 
194 
29 
158 
iv 
138 
199 
200 
ἘΠῚ 
eo 
205 


204 
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X44 DRAW THE MESH SURFACI OF THE GRAPH sw 


ERE E 555872 ο ولت‎ (٢٢ στ ασ η CEV O a Gers Ge 2a, Oe Sy IDIV, O, 


Do HARTI 2897 ت2‎ RU a MAS VERT SA) 

*+*+* ANNOTATION OF THE GRAPH ههبج«‎ 

ΕΗ m وژ لا 177 و 219 :لا وص لا‎ 9:1 τῷ) 

GO P ٢٢٧٧/٧٢٧٢ 22-90 99920 9 21 2221719 (ګک ,لا .لا‎ 

5: SNR Esa =; S=; a > u SIE Sah JN G) 

BALL RUMIEHER CS O IS O OE EEE ο ος 5) 

DM — 0201, 17730) * Ἐν ο 

BALL = Kilo Ho ο ο 2 Cl, 1) -DY, XR, YR) 

SALE 5هت‎ ٤ل‎ ۷۱١ لت )للا‎ fF ὁ, ὃς τὸ 

CACC 353۷1180 ο ο "بت ولا وا! لا‎ ## = ORIGIN? , 9.0, 10) 
ENDIF 


τ ΞΕΡΕΙ, : “ΞΕ. =m Co has S) : 
CALE SYMBOL كت بادا ات مرکا( (انرکا))‎ ε-ο DATA FIELD’ , 9.0, 14) 


**+*** QUTPUT THE GRAPH *#+#% 
CACC PLOT (9.0, 0.0.9399) 


WRITE (*, 416) 

READ (#*, 200) ANSWER 
IF  ((ANSWER .EQ. 'Y') 
WRITE (*, 417) 

READ (*,200) ANSWER 
IF ((ANSWER .EQ. 'Y^) 
STOP 


-ωπ. AANSWER «ΕΠ. ντο στο το 


ο. CAMELS τοι yy )) GOTG to 


ΞΟΒΜΑΤ/(Ω5ῷ) 
FORMAT (A) 


3 FORAT Z E0 OES tE BS KS رز هد‎ 


FORMAT () 

FORMAT (/, SX, AGO) 

πα οτι τν σι ση Sel)? 46x." (AZIMUTH 250.0)", 7) 
πο στ. SX CBI οσο Oo)? ASX, ° (AZIMUT 240.0)°, 7) 
FORMAT (10(F7. 2, 2%) ) 

FORMAT (9X, \) 

FORMAT (/, 5X, ! SEND GRAPH TO THE PRINTER(Y ar N) 
FORMAT(/,SX, 'DIMENSION QF OUTOU MATRIX(1 to 3 
FORMAT (/, 5X, ' ORDER OF TRANSFER FUNCTION(O to 4 
BROT SX BSI bis ee SD 
FORMAT (/, 5X,’ AZIMUTH(O.0 to 360.0 DEGREES): ὃς) 
FORMAT (/, 5%, °>ELEVATION(90.0 ts -90.0 DEGREES): °,\) 
FARMA EX τπιπο τα σος —Y$)$:$ ὃς 

FORMAT (/,5X,’2,4 OR 8 SUBGRIDS: ',N) 

FORMAT (/,5X,’ TITLE OF GRAPH(UP TO =O CHAR): °,\) 
FORMAT (/,5X,’DO YOU WANT TO CHANGE PARAMETERS 7? ',N) 
FORMAT (/,SX,’DO YOU WANT TO REPEAT THE PROCESS ? ',*) 
FORMAT (Z7, SX,°DO YOU WANT TO MAKE GRAPH 2 ',N) 
FORMAT(/,5X,’DO YOU WANT TO CORRECT ? *,\) 
FORMAT(/,SX,'DO YOU WANT CONTOUR MAP ? °,\) 
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AZIM 
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CTE xk 
عا اا‎ 
DY 


El SEMEN 


mE 
FLOAT 
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1400) 


IMGDEN 
IMGNOM 
IMGPSR 


IPROJ 
IT 
ITRIM 
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Re 

RLDEN 
RLNOM 


RLPART 


SIN 
Sü 
Sree 


VERTEX 


w1 
We 
AL 
XR 
XUPR 
Y 55m 
YR 
YUPR 
Z 
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Type 


REAL 
CHQamR*1 
RERL 
REAL 


CHAR*E0 
REAL 
REAL 
REAL 
REAL 


INTEGERw*Z 
INTEGERS 
INTEGER#2 
INTEGER*Z 
REAL 
INTEGER*Z2 
INTEGER*2 
INTEGER*Z 
INTEGER*z 
INTEGEAR*Zz 
INTEGER*2 
IM eSE RS ee 
Ay WSS ee 
INTEGER*z 
INTEGER*Z 
NESE Ra 
INTEGER*Z 
REAL 
REAL 
REAL 
REAL 
REAL 
REAL 


REAL 
REAL 
REAL 
REAL 
REAL 
REAL 
REAL 
REAL 
REAL 
REAL 
REAL 


18214 
IB 

18254 
18130 
13206 


18082 
LESIE 
181 76, 
18170 
J 3€ 3 بد پد‎ 
18074 
18078 
16e20 
18090 
13080 
18 152 
17850 
17952 
1515 
11850 
ALSS 
18076 
TIS 

534 

786 
1821 - 
LI 


3 3F 3F 3F 3F 


18124 
11786 
18128 
18140 
18054 
18248 
18062 
το τα 
στο 
18066 
1178 


INTRINSIC 


INTRINSIC 


1171 1115 8 
INTRINSIC 
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το -ν 
Ζι ον 
ΖΜΩΧ 
ΖΜΙΝ 


Name 


MAIN 

MESHS 
NUMBER 
PZD2D 
PLOT 

αι ΤΕ 
SYMBOL. 
4INDOW 
ZCNTUR 
σι Εν ΕΙ 


REAL 
REAL 
REAL 
REAL 


Type 


Pass One 


cO 
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n D 


شت د لل 
183070 
18115 


Laine 


Size 
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PROGRAM 

SUBROUTINE 
SUBROUTINE 
SUBROUTINE 
SUBROUTINE 


SUBROUTINE 


SUBROUTINE 
SUBROUTINE 
SUBROUTINE 
SUBROUTINE 


Errars Detected 
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FORTRANT7 


1207516 
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MEC 


DATA 


0-ه 











156 


|) 
αι 
. 99ص‎ 
اا‎ e 
- 
ri) 
ας 


ΕΠ Εν. 


πο η. 


[ N 


رل 


ΙΒ] 


سر 
B‏ 


-- 


tt: 


APPENDIX C 


Jage i 
)( 32-2 -ع‎ 3 
ers sce 
ο τοι I 7 EE SEE ARN ANA Vo, eU OSES 
1 SSTORAGE: Ξ 
2 $SPQGESIZE:Z8 
4 E dE JE JE JE Je JEJE خي خي‎ ie i JE JE پد‎ ib i د‎ ib i جد چېي ی ی ا‎ a 4 چو خد ځد چی چي چب‎ JE dE JE Jb De JE چې خي ځې سی خي‎ JE de te Ab ue d خي‎ P ځي خد سي‎ B te مک‎ 
Er * 
E * THE PURPOSE OF THIS PROGRAM IS TO COmPUTE AND GRAPH TRE - 
E + EQUATIONS OF ROBERT 2. ROESSER IN THE “DISCRETZ ΒΤΩΤΕΞΓΞΡΩΟΕ * 
SIE * MODEL FOR LINEAR [MAGE PROCESSING". IT TRANSFORMS ALSO THE + 
2) ® + OUTPUT MATRIX Y ACCORDING YO FOURIER ANALYSIS. * 
roO c * * 
LI C * EVANGELOS THEOFILOU * 
12 i Jb JE JÉ Je JE Je JE Je JP HIP JE HIP AP کب جد غد کب ا ج جیب جب‎ ABO جد جیب‎ ÀJ6 ÀP JP ÀE JC Je Je ند ځد چو ی کو ی جد ج‎ AE ځد خي‎ e ا‎ e 6ل ی‎ E ځچې‎ i ځد دل‎ ÁÀE 6-5 ΞΕ e e 6-6 خي خي‎ 
13 8 PROGRAM 2D-DATA-FIELD 
14 
= *3«333 YARTARBLE DECLARATIONS +--+ 
16 REAL R1 (26, 26), RZ (26, 25), S1 (25, 25), S2 (26, 25), 
2 * FRI (2), FRE(Z), TRM(4, 4), IV 04), OV (4), IMGPART 
15 CHARACTER*L ANSWER 
15 š 
5 κκ VARIABLE DECLARATIONS FOR PDLOTS8O8 s 
E CHARACTER*«zO CTEXT 
22 COMMON /WORK /Z (26, 26), ZF (26, 26), ZLEV(26),_LDIG(26), 
23 * LWGT (26) , MASK (5000), VERTEX (16) 
24 
c DATA ως. ο Avie ο ος XUBHR/O.TS/,YUPRX/7.0/, 
26 * ZLOW/ 1. OESS/, IPROJ/O/, NRNG/100/ 
E 
2 È p ie ie db cac d dde A db e HP π ο τν 2 RJ GRAM ee HE JP d sP 
So cC σος. OSH THE REQUIRED VALUES FOR THE mODZL *-—- 
S 10 WRITE (*, 4023) 
Se READ (τν ο) KK 
E k (Qa TOS DS TORS (KK Ol. 299) GOTO 10 
34 
E) DO 100 I 5 1,KK*1 
1 35 DE TOOT = 1, KKEL 
= my πι ιτ ο - θεὼ 
= 38 πο ο = 0.0 
2 So 59۱٩) ل٧ جح (ل‎ 0.0 
2 40 SE(I, J) = O. O 
E 41 100 CONTINUE 
42 
43 ο 191 تد 1ل‎ 1,4 
1 44 ΒΟ 101-7] - 1,4 
2 E πο ο. οὐ 
= 46 101 CONTINUE 
47 
48 στους in 
1 49 FS (To = 0.0 
1 ZO DUCI) سح‎ 0.0 
1 531 102 CONTINUE 


15; 


D Line= 


r^ 


Fr 


s= 
ss 
هه‎ 
ss 
sE 
s7 
-8 
z3 
&o 
61 
6.5 
65 
64 
65 
ΕΕ 
67 
68 
69 
7ο 
71 
72 
73 
74 
75 
76 
v 
78 
79 
ao 
81 
83 
a3 
84 
as 
ae 
a7 
aa 
a3 
30 
31 
92 
93 
94 
95 
96 
97 
38 
33 
100 
107 
102 


103 


104 


106 
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WRITE (* 
DO 


1 
EOS: 1 
WRITE 
READ 
CONTINUE 


WRITE (*, 
DO 104 I 
WRITE 
READ 
CONT INUE 


WRITE (*, 

DO LOS ٢ 
WRITE 
READ 

CONTINUE 


WRITE 
DO 106 I 
WRITE 
READ 
CONTINUE 


WRITE 
Εν ο ο = 1 
WRITE 
READ 
WRITE 
READ 


Ge, 
) 
ο 


WRITS (#, 
TRM(1,2) 
TRM (4, 1) 
WRITE (*, 
RERD (*, 
WRITE 
READ 
WRITE 
READ (3. 
TRMCIS 2 
WRITE 
RERD ) 
TRM(1, 4) 
WRITE 
READ (*, 
TRM Ce. So) 
WRITE 
READ (3. 
TRM(Z, 4) 


(*, 


711) ENTER 


211) 


err) 


*) 


409) 


*) 


1 


eo of 


$09) 


*) 


*) 


*) 


*) 


*) 


*) 


1 


(3, 604) 
(*, *) 


1 


(σε, 404) 
(*, *) 


1 


(*, 4035) 
(*, *) 


CE ND 


1 


(#, 405) 
QR) 


(*, 40) 


1 
1 
1 


(3:90) 


(*, 409) 


, KK 


' ENTER 
,KH 


' ENTER 
, KK 


' ENTER 
, KK 


(#, 211) "ENTER 


OG ies 
OV(3) 
aw iss 
OV (4) 


) TENTER 


αι. 
TRMCIIS 
* acus 
TRM(Z, 1) 
Y aj ol ¢ 
TEMP 


12117 
TEMP 


"Bets 
TEMP 


xe: 
Temp 


Micrasoft FORTRANTT7 


INITIAL CONDITIONS 


RIC TT 
ΠΊΝΕΙ, τ} 


INITIAC CONDITIONS 


πι. 1.575.” 
R2(1,1) 


INITIAL CONDITIONS 


δι τὸ ο 
SIMI. 


INITIAL CONDITIONS FOR VERTICAL 


"ester e MS 
Sa (1, 1) 


VALUES FOR THE OUTPUT VECTIOR(CH. 4) 


ELEMENTS OF- THE TRANSITION MATRIX (4. +) 


TEMP = ον ο ο iQ) 
(*, 409) 


TEMP + OV(4)*TRM(1, 1) 
(*, 409) 


TEMO + DV (S) *TRM (Z, 1) 
(*, 409) 


TEMP + OV(4) ¥TRM(2, 1) 


138 


SOR VERN EAC SI CR) 


Se (#. #) 


FC də TZ sS so Was ΙΙ; 


POR ποστ Ντο RSCHOSO" 


9 


9 


1 


fo fu fg fu [o fo fo [o e 


fis fu re 


fu fu 


b 


DJ 
n) 


عم هم وم مم عم 
(JJ (J) (jJ CJ CJ (J) O) C‏ 


Qo - O (n P CJ 


C^ ror 
+ 
سا‎ 


- 


t 
+ 
t 


142 
143 
144 
145 
146 
147 
148 
149 
S0 
151 

E 


سب 
-- 


κ. 


107 


198 


109 


110 


~ace 


2 
— 


„re - = ig 
Dg roS 


= ana = q a 
7 MIEresarft FORTRANT7T Vo CEOS AA 
τοις ο “σου i 
TRM(4,4) = DV (6) 


WRITE (s, 
ντο لل‎ 
WRITE €€,402) 'boo: : 

READ (æ+, +) IV(4) 

Απ ΙτΕ ἑ πώς) “bigs ” 

READ (*, *) TEMP 

IV(1) - TEMP + IV(4) *#TRMC1, 1) 
IV(9) * IV(A4) *TRM(CS, 1) 


LAA 20 VALUS PUN raa TNE VEE AR CHE RD 


[a 


a 170 
DO 207 1 1K 
DO 107 J = 1,KK 
Rit Dil SRI PD τετ لل‎ + TAMIL. ο κδι(τ + 


ΤΗΗ (194) «551. 0) + IV(1) κ} 
RÆ(I+1, J) 2 TRMCE, 1) *RI CI, J) 9 TRM(CO, 2)9*S1(CI,J9 + 
Pama. +7 ©5201, J) & TV C2) £0 
SI OT OSI SY 
S2(I,J*1) = RICCI, J) + QOV(3) *SICI, J) +* OV (4) #SZ(I, J) + ἵν (2) ει) 
U = 0.0 
CONTINUE 


WRITE (#, 205 
WRITE (*, 3 


) Txeee# INPUT VECTOR κκ) 
200) CIVCID,I = 1,4) 


WRITE (æ, ZOS) "s QUTPUT VECTOR oe? 
Werle ος ο ως «ο ος A) 


WRITE (#, 205) ' eee TRANSITION MATRIX “4448? 
DG ος 1 - هول‎ 

WRITE (+, 200) (CTRMCI,J),J = 1,4) 

WRITE (#, 210) 
CONTINUE 


eee FILL O's THE TWO DIMENTIONAL GRID OF CONTROL POINTS «*«% 
DO 109 I = 5 
DO 109 J = 1,26 
ZITAT) aon) 
CONTINUE 


DO 1160 I = 1,KK 
DO 110 J = پا‎ 
Z(1,J) = R1¢I,J) + OV(S)#S1 (1,5) + OV(4) *G2(1, J) 
CONTINUE 


WRITE (#, 205) '«eeee RL MATRIX ',KM,' X ",HK,? sex? 
pu a πα 

WRITE (#*#,300) (R1¢(I,J), J = 1,KK) 

WRITE (#*, 210) 
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D Linett 
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هم 


154 
133 
126 
1-5 
128 
tos 
160 
1Επ 
162 
163 
164 
165 
166 
187 
168 
l5 
170 
131 
172 
lm 
174 C 
19 
176 
157 
178 
1 د‎ 
180 
181 
182 
185 
184 
183 
186 
12۱77 ت)‎ 
138 
197 
{51 
TM 
132 
Las 
134 
132 
139 
137 
138 
133 
ZOO 
eol 
2032 
E E 
204 


ID 


Liz 


113 


115 


ni M2 IE SE Fl 


νο) 


PARAME TERS 


t 


3808 
ΞΕ 
Sees Se 


ΞΓΌΡΤΠΡΝΖΖ V3. 20 O 


JJ 3 396 9€ 3 


O 
x μα) 


ا د جد چڊ چ 


X 3 3€ 26 46 ° 


X ".KA," 


385 3€ 36 3€ 96 ? 


X EH s 


Go 21 


E E d دا‎ 7 


CONTINUE 
WRITE ο soOS) *xxcx RES MATRIX TRE, 
DE x le 7 — پا رل‎ 
WRITE (+, 300) (Retl, J), J = 144%) 
WRITE (#, 2109) 
CONTINUE 
WRITE (*,205) "'xwx 51 M ATR κ ο... 
ρο 1135 I - 1I. M^ 
WRITE (*,.500) πι κι οι UL ww TE 
WRITE (#, 210) 
CONTINUE 
WRITE (#, 2035) ‘x*eee SZ MA TR IX °K, 
DO 114 I = 1,KK 
WRITSNO*, 200) 3452 (15), J. = απ) 
WRITE (#, 210) 
CONTINUE 
κκ OUTPUT THE Y MATRIX s 
WRITZ (X. 202) ΜΕΣΑ Π.Ι πα EID 
WRITE (*, 212) 
DOE 1157 Ἱ ο... 
WRITE (4, 00) CZ CTS ο πο 
WRITE (κι 210) 
CONTINUE 
White Cee ts) 
WRITE (*, 419) 
READ (#, 200) ANSWER 
IF ( (ANSWER «ΝΕ. ?Y?) J AND. (ANSWER . NE. 
33 ASK THE PARAMETERS FOR THE GRAPH wTw**** 
WRITE (#, 210) 
μπες te.) "cee E NT EDT 
WRITE (#*, 410) 
READ (#, #) AZIM 
WRITE (*, 411) 
READ 5) ELEV 
WRITE (#, 413) 
READ (#, *) ITRIM 
WRITE (*, 414) 
READ (*, *) τοι 
WRITE (#*, 415) 
READ (9:133) EXT 
WRITE (*, 451) 
READ (3, 200) ANSWER 
3*3*3*3* INITIALIZE PLOTS88 xxx» 
IF ((ANGWER .EQ. *Y') للا‎ ΙΡ επ Ευ. τ 
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ED 
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GO 
208 
223 
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211 
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213 
214 
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216 
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218 
219 
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fo fo [O fe fo fu fs fo fü 


LUE 


38968 


03-26-52 
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کت کب و( اه 
Micrasoft 20817۲۹۸٩7 7 Vi. dU 02/84‏ 7 
CALL PLOTS(0, 0, z)‏ 
So Se‏ 
CELLS Peers ((, 22, 32)‏ 
ENDIF‏ 


CALL WINDOW (XLOL, YLOL, XUPR, YUPR) 


**4 DRAW THE MESH SURFACE OF THE GRAPH x 
CALL MESHS(Z, 26, 26, KK, KK, AZIM, ELEV, 0.5, 0.5, 9. 25, 6.5, IDIV,9, 


* 5, TPROJ, 1, ZLOW, 3, ITRIM, MASK, VERTEX) 


eeeee ANNOTATION OF THE GRAPH teres 

εσυ ο Παν D. So Z=. σαι Ετεκτ ο. ο, 5ο) 

ο SMES τν ο τα eee ce -e DATA FLELCD*..0.0, 14) 
ΕΕ ΙΕ ο ο ο ο ο ο  Ασ  ΜΙΤΗ: ۰ل وي‎ 9 του 
ΕΕ ο ο οσο. ος ο, δις 2هو‎ Mo Oo G 2) 

Γρ ιο πο tas. OO. Ome. 1 ع‎ ντι ,و‎ 9:0. 10) 
COLL NUMBER ( 33370..929..0, 0. 2, GLEV, 0.0, 2) s 
DY = (Z(1, 1)/3O.O) æ ELEV 

στ ο «ου κα να)‏ 9 دل 0 ل2020 )نې 

Επι. 2۹501 ας YR OES, 0.0, 1) 

Cakes SY Mego, Osc. ate = ORIGIN’ ο ο 1ο 


**+*+ OUTPUT THE GRAPH sx 

EBLE BEAN CO OO O 999) 

WRITE (5.416) 

READ (#, 200) ANSWER 

ie CCRNSGER اتان‎ ο ο (δα (EO. ΤΥ} GOTO 20 


WRITE (#*, 418) 
READ (*#, 200) ANSWER 
Lm ( (ANSWER EQ. 'Y') .QR. (ANSWER .EQ. *y’)) THEN 


wc) FILL 0's THE TWO DPIMENTIGNAL GRID GF CONTROL POINTS xx 
ρα πιο ۸ 2 تح ولا‎ 
DO 116 J = 1,26 
στ ο) =O OG 
CONTINUE 


ΖΕΜΩΧ - -5.9Ε50 

ZFMIN = 9. 9E20 

DK = (KK - 1) / 2.0 

8 = 3.141532 

DO 117 M = 1,KK 

DO 117 N = 1,HK 

RLPART Ó. Ó 
IMGPART = 0.0 
DO 118 L = 1,KK 


DO 118 ٤ 1, KK 
۳٣ ٤ ος EAS E EE © (l= ines M—PK—1) 7K) 
ΕΠ1 ΤΕ) Αι ακοτ ος ο ετήσιο KK) 
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fu fo ٢ Fj fj [Û û Û $È FF دود د د د‎ Ü 


Lines 
256 
c 
258 


Ao 
= 


ceo 
ο] 
oO 
ce 
264 
cog 
eat 
267 
268 
ΞΕ 
cy 
art 
eye 
evs 
274 
279 
= = 
ei? 
=o 
ne, 
28o 
281 
cT 
BS 
284 
235 
ES 
SOV 
288 
289 
290 
eas 
ese 
eem 
294 
ig 
σα 
E 7 
25e 
S 
a OD 
301 
302 
303 
2504 
το 


306 


115 


117 


113 


ت .ت۷ 7 2۱۳۱۹۳۸7 mMicmrasoft‏ 


με. COS (-2*P* (K- 1) € (N-DR-1) /KK) 

PRECES SIN (SSO E CES LD E CNS DE DAES) 

RLPART = RLPART + Z(L, K) #(FRI (1) #FRS3 (1) 
ERI CE) 2م حم‎ EDD 

IMGPART = IMGPART + I(L; KI#F(ER IG) sp a= (S) 
E RI COD c= = 1 D) 


CONTINUE 
ZF (M,N) = SORT (RLPART#*2 + IMGPART **2) 


D 


( 


ZF(M,N) .GT. ZFMAX) THEN 


ZFMAX = ZF (M, N) 
ENDIF 


IF 


(ZF (M, N) SDT ZEMIN ΠΕΝ 


ZFMIN = ZF (M,N) 
ENDIF 
CONTINUE 


ΕΕ OUTPUT THE ZF MATRIX w**** 


WRITE 
WRITE 


Cx, 


(ας 


νι ΝΕ 
WRITE 
WRITE 
CONTINUE 


WRITE 


(3* 


205) +++ FOURIER TRANSFORMATION  ',HK 
,212) 

= 1,KK 

(%, 300) (ZECI; J), Ta KK) 

(#, 210) 


913) 


WRITE (*, 419) 


READ 


Dm ( (ANSWER ΝΕ. 'Y') .&AGND. (ANSWER .N&. *y’)) GO 


Á*3* 3t 363€ 
WRITE 
WRITE 
WRITE 
READ 
WRITE 
READ 
WRITE 
READ 
WRITE 
READ 
WRITE 
READ 
WRITE 
READ 


3e3* 3€ 3e 3€ 


(*, 


as 
(* 


IN 


200) ANSWER 


K THE PARAMETERS FOR THE GRAPH κε 
210) 
3 


που νου ἙἝ μοι ο πο ο. OS p P J Rea (= 
; ^10) 

١٢ AZIM‏ و 

5 411) 

, *) ELEV 
3511) 

; x9 ITRIM 

3 414) 

IDIV‏ )* ې 

> 4135) 

31395 ΕΑ 
» 451) 


, 200) ANSWER 


ITIGLIZE PLOT88 xx 


Pan 


.. 


ἡ 
ft fu m 


-. ο 


τη. (Ὁ 


a 
BL 
mC) 
O 


Ε [i Û iN 


ο 


κ μμ” ومر‎ 


TORES 


T ER Sas s, 


IF «(ANSWER .EQ. ¥) .OR. (ANSWER .EQ. 'y')) THEN 
επι 660756 O) O) 


ELSE 


CALL PLOTS 99 93) 


ENDIF 
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D cines 
S307 
308 
e 
3 193 
ult 
Se 
ld 
314 
Bis 
216 
1 17 
5 
1 E 
E 
ul 
See 
G = 
324 
22S 
Se 
E 
Ges 
pers, 
339 
Sal 
Soe 


— 
اه ا ا 


234^ 
22 
coo 
7 
338 
g, 
240 
341 
342 
343 
344 
245 
346 
347 
28 
243 
uz 
ο 
ات زی‎ 
Su 
E 

SE 
Sa 


σα Ἐπ 
w J 


fo 
fo 


136 


139 
5ου 
ΞΟ5 
210 
e 


lu 


WRITE (#, 420) 


READ (#, 200) ANSWER 


πια σας 


CALL WINDOW (XLOL, YLOL, XUPR, YUPR) 


IF ( (ANSWER . EQ. 


τω 


OR 


. (ANSWER 


DLEV = (ZFMAX-ZFMIN) /FLOAT (KK) 
αρ σα ο. ο ο ας δε σεν κ.) 
DO 136 I = 1,KK+1 


LDTOWVI) = 2 
LWGT (I) 1 
CONT I NUE 


«EG. 


tO) 


PORTRANZ 7 


THEN 


- { } 


la fa 


~ 
i è 


> D DY 


- 


] 7 


fu 0 


fu [u Ch 


Lr 


.- 
٩ 
-— 


CALL ZCNTUR(ZF, 26, 26, KK, KK, 0.5, 0. 5, 9. 25,6. 5, ZLEV, LDIG, LWGT, 
CALL SYMBOL (S.5,0.0,0. 2, ' CONTOUR MAP',0.0,11) 


ELSE 


«33 DRAW THE MESH SURFACE OF THE GRAPH sx 
CALL MESHS(ZF,26,z6,KK,MM,AZIM,ELEV,0.5,0. $5, 8. 25, 6. E, IDIV, O, 
3, IPROJ, 1, ZLOW, 3, ITRIM, MASK, VERTEX) 


*+*+++ ANNOTATION OF THE GRAPH «ελ. 
E AEE ΕΓ CES وک ولا رک‎ i". mZIMUTP: 
caL MUER EEE O) EEE ο ο ο, AZIM, oO. UV, 2) 
ο ٨۴ 5٢ 7 ου ο ος Ε͵ι ενατιΟοΝ το. ὃ, 1) 


GREEWRUMBER CJa9. 0, 2232. 0, 


DM 2) ولا‎ ον πω. ο) 
SALE 22020 1.9 لل‎ Z= ντου απ Υπ) 
EIB S YE (E) PIE ولا :۱9 ولس وحصكت ولا‎ 1) 
ΓΗ: ΜΕ ο ο I. Ou 


ENDIF 


κ ELEV 


ο 0) 


O. 2, ELEV, 0.0, =) 


na ORIGIN? OT OTIO) 


EA EE SYMBDOLUISUSS: 75. 0. 25S, CTEXT, 0. 0, 
Coll sy MEGLC KG eee-a, Dee, '2=—D DET*, Oo. 0, 7) 


κκ OUTPUT THE GRAPH κ.ε 
τα αι ο ο ο ου 53 9) 


WRITE (#,416) 
READ 


(*#, 200) ANSWER 


IF (CANSWER .EQ. *Y’ 


ENDIF 
WRITE (#, 417) 
READ (#, 200) ANSWER 


IF ((ANSWER .EQ. 
STOP 


FORMAT (ASO) 
FORMAT (A) 


π᾿ 


m 


ο OR. 


R. (ANSWE 


(ANSWER 


FORMAT (/, 18X, AZ9, 12, AS, 12, AS, /) 


FORMAT () 
FORMAT (/, SX, ASE) 


R 


= 


. cC 


eo) 


EU. 


D. 


384 12)8) 


*v')) GOTO 


GOTO 10 


FORMAT (/, 2X,’ (AZIMUTH 320.0)',46X,' (AZIMUTH 220.0)", /) 
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20 


N هوه‎ | 
© (js Ὁ 


ro ch st 


í 


D Exnes 


Nam 


= 


ον 

Soo 
S60 
Got 
OS 
SES 
564 
265 
266 
3267 
368 
263 
270 
S71 
ym 
37 
374 
= 
376 
377 
378 
oro 


ANSWER CHAR*1 


AZIM 


CGS 


CTE 


DK 


DESY 


DY 
ELS 


V 


FLOAT 


PRI 
ΕΠΕ 
1 

τοι 


ν 


IMGPAR 
IPROJ 
ITRIM 


IV 
J 
K 
KK 
L 


HOE 


LWG 


T 


1 7 Micrasart FORTRQN77 
213 FORMAT (/, 2X, 7 (AZIMUTH 65050)*- 46X AZ TT 7 
290 FORMAT C1O(F7. 2, 1X) ) 

400 FORMAT (9X, \) 
403 FORMAT (/, SX, "DIMENSIGN GF OUTPUT (A—ltca eae 
404 FORMAT (SX, AS, Ie, AS, N) 
4035 FORMAT (SX, AS, IŽ, AJ, \) 
406 FORMAT (SX, A, 12, A4, \) 
407 FORMAT (SX, Aa, I2, IZ, A3, \) 
408 FORMAT (SX, AS, le, AS, N) 
409 FORMAT (SX, AS, \) 
410 FORMAT (7, Sx, °AZIMUTH(O. 90 τε 260.0 DEGREES) 2) 
411 FORMAT (7, Sx, * GEEVel LON C20. Oo صن‎ πο ο Εβρ ΕΕ ο... ss) 
415 FORMAT (/, SX,’ TRIMCOSNO, 1=Xs, 2=Ys): 7, N\) 
414 FORMAT (/, SX,'°2, 4 GR 8 SUBGRIDS: 'S X9 
415 FPORMAT(/, SX, TITLE GF 5:929) 5 το ο ας ο κ 
416 FORMAT(/,S9SX,"DO YOU WANT TO CHANGE PARAMETERS” 7,\) 
417 FORMAT(/, SxX,° DG YOU WRNT TQ REPERT THE PROCESS”? ',N) 
&#18 FORmAT(/,SX, *'DD YOU WANT FOURIER TRANSFORMATION ? ',N) 
4193 FORMAT(/, 5X, *'DD YOU WANT TO MAKE GRAPH 7? ',N) 
420 FORMAT(/,5X,° DO YOU WANT CONTOUR MAP ? 3 ,N) 
451 FORMAT (/, SX, ° SEND GRAPH TO THE PRINTERC(CY or N): ',w) 
END 
Type Offset P Class 
11060 
RERL 1 
INTRINSIC 
CHAR*#20 11074 
REAL 11114 
REAL Liles 
REAL 11994 
REAL 11066 
INTRINSIC 
REAL 10882 
REAL 19830 
INTEGER EE 1099€ 
INTEGER#Z IBNOE 
REAL 11142 
INTEGERS 1095€ 
INTEDER*Z 11070 
REAL 19898 
INTEGER*2 10364 
INTEGERS 11154 
INTEGER*2 10924 
INTEGER*Z 11146 
INTEGER*e2 Sl /WORK 7 
INTEGER*Z 35364 /WORK  / 
IMTEOGER*S 1112: 
INTEGER*:-Z 3616 /WORK  / 
INTEGER*G 1119 
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Pace 8 
υπ ef ea 


20:42:32 


VS. E0 


2 / 84 


D Lines : 7 

NRNG INTEGIÉD*c 1093952 

ον RESL 103124 

p RERL 1111Ə 

R1 REAL 3 

Re REAL ΞΕ 

RLPART REAL τι 58 

St REAL 2410 

Sc RERL 2114 

SEN INTRINSIC 

SQRT INTRINSIC 

TEMP REAL 10226 

TRM RERL 10818 

UJ REAL 11900 

VERTEX REAL 11616 /WORHK  / 

XL OL RERL 10329 

XR REAL 11098 

XUPR REAL 10558 

YEO REAL 10324 

YR RERL 11102 

YUPR REAL 10942 

2 REAL Q /WORK ٢ 

zF REAL 2704 /WORK  / 

ZFMAX REAL 11155 

ZEMIN REAL 1 11 

ZEEV REAL 5408 {WORK / 

ZLOW REAL 10946 

Name Type Size Class 

MAIN PROGRAM 

ME SHS SUBROUTINE 

NUMBER SUBROUTINE 

P = = SUBROUTINE 

PLOT : SUBROUTINE 

PLOTS SUBROUTINE 

SYMBCL SUBROUTINE 

WINDOW SUBROUTINE 

WORK 11680 COMMON 

ZCNTUR SUBROUTINE 
SUBROUTINE 


LEVEL 


Pass One 


No Errors Detected 
3723 Source Lines 
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MASS =i 


U 
"apos 
| wo 
li. U 


(D (J UD 


-- 
ME 
ل‎ —— 


fi (ñ u) 


- M 
٩ 

— ~ 
- 

-— 


art 
F 


fo fie qa 


ci 

1١ 
Ό 

له 

4 

al 

Db 
q 
" 
< 
(J 
[a Du 
`, هوه‎ 


Oo Lines 


GJ (J GJ ω fu هم‎ 


fu 


tr 
o O O YN gi 5 for 


هم هم مهم هم 
(J D) —‏ + 


مې 
(η‏ 


هم مم مم مم 
u) On Mm‏ 


4 


[J fi fà PO fa fu fJ [o [0 fi 
) Í) — 


Uu) (€ cO £0 PC 


40 


42 
43 
δῷ 
45 
46 
47 
48 
49 
ΞΟ 


= 
اب 


APPENDIX D 


Microsof? 


Dage 
03-26-5839 


٩ 

- 
چم 
—- 


19:32:06 


FORT GONYV ود کر‎ 


02/94 


T-t خد نې خي خي خي‎ ie 3E JE JE JP JR JE i ie a يف‎ HF e <b خد‎ NË JE HO JP EP HF HF ae i i et i tt i i خي‎ ie ie E NP. خد‎ E e ie i i -ᾱ- -ᾱ- 


TO 
ΠΕ 


COMPUTE AND GRAPH THE 
“DISCRETE يب 5۴ - ع5185‎ 


TT TRANSFORMS ALSO THE 


EVRHNGELOS ΠΠΕΕΠΙΕΞΗ 


ΕΕ‏ اد اد جد جد سل د د د د د جد جد جد جد يي خي خي اد اد جد جد جد اد اد جد اد اد اد د ج جد اد جد جد جد اد د اد جد د جد ΞΕ ΠΕ‏ جد اد ae JEJE JE Je JE EP ie i‏ چيا نه سي د خي يه سی یی ځا 


RL CED RE CED PRINS te) 1 ) ١ ۰( ره و‎ ٢۳و‎ 7115277 


- 


/WORK /2(26, 26), Ze (267 2607 202) (26) DIS (Se) 


ALOL/O.O/7, YLGL/O. O7 ,AUPR/S. S74 μα... οὐ, 


JE 3t Je JC Je JC te i ie ie ib ie سي‎ 


1 z 
SBSTORAGE: = 
SPAGESIZE:58 
نا‎ 
[z * 
E * ΤΗΕ PURPOSE OF THIS PROGRAM IS 
G * EQUATIONS OF ROBERT 2.ROESSER IN 
ë * MODEL FOR LINEAR IMAGE PROCESSING”. 
Ç * OUTPUT MATRIX Y ACCORDING TO FOURIER ANALYSIS. 
C ae 
[s * 
C 
(5 PROGRAM 2D-DATA-FIELD 
C κκ YORTABLE DECLARATIONS #++#+ 
REAL R(26, 28, 4), 31 (25, 28,4), 52 (26, 28.4), 
* 
CHARACTER*1 ANSWER 
C **+++ VARIABLE DECLARATIONS FOR PLOTSS wee 
CHORACTER*ao CTEXT 
COMMON 
* LWGT (28) , MASK (2000), VERTEX (16) 
DATA 
* ZLOW/ 1. OESS/, TPROJ707, NRNG? LOO7 
= T چب‎ eb ie -ᾱ- -ᾱ- Ib AR د يا څد‎ Mm A IN ο RQR O G 5 ٥ 
0 


19 


[9 


ω 


100 


εκ. ASK THE REDUISED VALUES FOR THE MODEL 


WRITE (#, 401) 
READ Qe) N 
IF COMIS INDE 0) 
WRITE (», 402) 
READ (*, *) m 
LIF «Πτι. 
WRITE (*, 403) 
READ (*,*) KK 
IF (KK .GT. 


«08. (Ν 
. OR, (m 


Em 529019 د‎ 
DO 100 I = 1,HK-*1 
DE IOD s لح ول‎ 
DO 100 L = 1 
που το πο 
τι CO. 
Se el poi E= 3 
CONTINUE 


DO 101 I = 1,N+2*™ 
DO 101 J = 1,N-*cZwM 
TRMCI. JO) - 0.0 


. GT. 


ST. 
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& ) ) 


( زې 


GOTO 10 


57 


* 3F 


f £ *‏ £ ې ې ې 


[ [Ὁ [ὺ -- fu fu fu [ὺ [ὐ [9 -- 


هم هم هم هم 


مهم هم هم هم 


مې 


E t هدم‎ 


920 
231 


= 
= 


34 


ET 
-ὄ 


36 
2 
e 
a2 
100 
101 
102 


101 


1o2 


155 


104 


LOS 


106 


107 


108 


103 


- 
CONTINUE 
DO IOS TIT SS INEZ CN 
τν) = 0.0 
OV (I) = O. O 
CONTINUE 
WRITE Cog ELLI ENTER INITIAL 
DO 105 T <> YOKE 
DO 103 J = 1,N 
WRITE ος 4O4) απ IS, (1, ? 
READ (x (خ‎ ٧) D. Τ} 
CONTINUE 
WRITE τ": ۷۱5٧/1257 01 1 11 
DO 104 I = 1,KK 
DO 106 J "ولتك‎ 
WRITE (#, 405) 
READ ο ο ol Ol. 15) 
CONTINUE 
WRITE جح‎ 211 RENNER INITIAL 
DO 105 I = 1,KK 
DO 105 J = 145m 
WRITE (#, 405) ٤ 
RE AD ον ο Sa Cl, (ل‎ 
CONTINUE 
WR Pewee) ° SNTER VALUES 
νει .لاج‎ 
ο τος I1, 
WRITE (*,408) a(0!,I,?): ) 
READ (FE © TVNET) 
TRM(N+I,N+1) = σπιν (Ντ) 
CONTINUE 
DOO تن از‎ 1,.™ 
WRITE Ca GOOD BDC ول‎ (۸:7 
READ (#, #) IVC(N+M+I) 
TRM(N+M+I, N+1) = -IV(N+M+T) 
CONTINUE 
Ρο 109 1 لولح‎ 
Π.Ε 3 1.0 
CONTINUE 
ΞΘ ΤΠ 
TRM (N*#M+I,N#M+I41) = 1.0 
CONTINUE 
KRITE Oe ENTER 
DO 110 I = 1,N 
EHRIIECCOE, 405) aC’ Siew) = + 
READ (#,*) TEMP 
TRM(1, I) = -TEMP 


FOR 


Pane e 
09-26-85 
19:32:06 


Microsoft ΕΟΠΤΠΙΩΝ V2.20 02/84 


CONDITIONS FOR HORIZONTAL R(#. #)? 


ele ο, 


CONDITIONS FOR VERTICAL SIRE. RD ' 


oN OCT IIS 


FOR VERTICAL S2(#.#) = 


THE INPUT VECTOR (t. 8) 1 


ELEMENTS OF THE TRANSITION MATRIX (#. 8) : 
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0 


— 


fo fo fu fe 


fu fu fu fue 


هم هم يم عم 


t + t O (J (J (Ü Q T +f T f j T +j G fe 


Line# 
DO 
104 
1035 
106 
107 
108 
103 
115 
Iq 
1 
115 
114 
LIS 
ΤΕ 
1 ٣ 
118 
112 
1 9 

c 


em) 
cc 


123 
124 
125 
126 
pe? 
128 


130 
Tl 

ae 
TSS 
134 
155 
135 
ر‎ 
τα 
139 
140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
L6 
151 
1 
1-5 


Jage 


DOST ESAS 
12: 32:5 
7 microsoft FORTRAN77 V3. Z0 0/84 
TIO CONTINUE 
TRMC1,N#1) = =1. 0 
DO 111 I = ابت‎ 
TRMCI.I-—-1) = 1.0 


111 CONTINUE 
DD 112 I = 1,m 
DO 1124099 
WRITE (*,407) "aC, J, I,? 03 7 
READ (*,*) TEMP1 
TRM(I*-N, J) = TEMP1 + TRM(1,J) * IV(N+I) 
CONTINUE 
DO 112 I = 1,M 
pO 115 7 = ME 
WRITE (8,407) DO TOTO 
READ  (*,*) TEMPI 
TRMCI+N+m, J) = TEMP1 + TRM(1,J) * IV(N+M+I) 
113 CONTINUE 


(0 


11 


WRITE (+, Z211) °’ ENTER VALUES FOR THE OUTPUT VECTOR(#. #) 
WRITE (*, 409) *b(00): " 
READ (*,*) TEMP 
OV(N+1) = -TEMP 
DV (N+m+1) = 1.0 
DO 114. 1[7= a8 
WRITE (FE GOS ες, τ, os 
READ (#, ¥) TEMP1 
OV(I) = TEMP1 + TRM(1, I) * TEMP 
114 CONTINUE 


)را 
DO 115 I = 1,KK‏ 
DOM ISEE epee‏ 
DO 116 II - 1,N-Z2«m‏ 


IF (II .LE. N) THEN 
DO 117 JJ $9 1,١ 
ΙΕ Ο.Ε. πιει 
R(I+1,J,II)=Rt(I*1, J. ID Παπ; Τη 
ENDIF 
IF ((JJ .GT. N) .AND. (JJ .LE. N+m)) THEN 
R(I+#1, J, 11) sR(I41, 35, II) -TRM(II, JJ) *S1 CI, J, JJ-ND 
ENDIF 
117 CONTINUE 
R(I+1,J, II) =R¢I+1, J, II) + IVCII) * U 
ENDIF 


IF («II .GT. N) .AND. (II .LE. N+M)) THEN 
DO 118 JJ 5 1,N-*Z*M 
IF (JJ .LE. N) THEN 
Si(I, J+#1, II-N) = S1i(I,J+1, I1I-N) + TRM(II, JJ) . 
* Pett, Jd) 
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لال + ++ r‏ خڅ MC‏ ())() خھ خښ خ خڅ خڅ خ خم خم خم 0 fu fut) C) ) E.‏ 


fu fù 


fer 


Lines 
s 


135 


156 
L57 
158 
159 
150 
161 
162 
163 
164 
165 
166 
167 
168 
169 
170 
171 
iS 
NOS 
174 
173 
176 
I 
178 
179 
180 
181 
185 
185 
184 
185 
186 
187 
188 
189 
190 
191 
192 
193 
194 
195 
196 
197 
198 
199 
ZOO 
201 
sog 
203 


204 


115 


GIS 


Sage + 


Q32-z5-az 
πο. 5 
7 f Microsoft FORTRAN? 7 Y3. 20 0/8- 
ENDIF 
DESIT رلا‎ 2۱9 (JJ «LE. N*M) ) THEN 


SL(I,J+1,II-N) = S1(I, J+1, II-N) 


+ TRM(II, JJ) * 


ιτ... J ND 


ENDIF 
CONTINUE 
τ: لل‎ -)٤( لح‎ ل٢‎ tl, LI-N) + ΙΙΙ) ευ 
ENDIF 
IF (II .GT. N+M) THEN 
88 لص ولګ ٢ل دحل‎ 
TE OIIOLE. N) THEN 
ت‎ 7111 ٧٢٢٢ ο S του οι τπτ 5: ΙΠΜΕΙΤ, Τη) 
κ ποτ. τα. 
ENDIE 
IF ETOT NNR ANDET ~.LE. Nem) THEN 
Zka Jilg = SS TSI ETN MD ο R Gr FOB To 
* Si¢(I, J, JJ-N) 
ENDIF 
IF (JJ .GT. N+M) THEN 
Soe ια = Sel, orl, TIK-N=M) + FRMCII, JJ) 
E GEL. J J N= 
ENDIF 
CONTINUE 
SEa Jelg LSNA IS Ισ ση و‎ 1٨ ۴1 «* U 
ΕΝΟΙΕ 
ONTINUE 
Quz 0.1 
CONTINUE 
WR τα το) Έντο INPUT VECTOR 33333" 


WRITE (*#, 500) (IV(I),I = 1,N+e*m) 


WRITE (#, 20S) "sex QUTDPLIT VECTOR sex? 
WRITE (*,Z00) (OV(I),I = 1,N+2*m) 
WRITE (#, 205) ex TRANSITION MATRIX sex! 
DO 120 I = 1,N-*2*M 

WRITE (=, G00) (TRMCI,J),J = 1,N+2=«m) 

WRITE (*, 210) 
CONTINUE 


κκ FILL O's THE TWO DIMENTIONAL GRID OF CONTROL POINTS #4 


DO 121 I = 1,26 
DO 121 J = 1,26 
LOI do 
CONTINUE 
DO 122 1 2 1,KK 
DO 122 J = 1,KK 
DO 123 LL = 1,N+2*M 
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i Û) (J (J CJ) CJ CJ (J (J (J) (C) O 


هم هم 
Γι) Γι‏ 
Dà) (J‏ 


0 


mee [u [à — 
LJ) [à + 

نا 
F Ln‏ 


[u 


Dy [O P [O D) PY) [ú D) [ó [ó 


.د 


Γι 


mee fy TO : 


خم 
m‏ 


n 


J (J CJ CJ CJ CJ 
n 5 (J b) لح ظا نا ن) حم‎ m (n + 


[g [O Da DJ [ú [ó Du fu Fo Fo DD F9 pà Fo [o fy [Oó 
C 
C 


[Γι 


me fo fue 


fia 


(n cn cn cn cn cn 


fo [fo [2 [o p fo 
(m £ C) [üu C 


(;΄3- 
12 
Z Microsoft -*-ORTRON?7 M3. 30 
Im s N) ΤΠΕΝ / 
تت (ل ول )2ه‎ σι τος + OV CEE eek eee ey 
ENDIF 


IF COLL GT. N). AND. (CORTE Nee 

C CT 1 ZCOI.32 o QVC) + ose ee 
ENDIF 
IF «DLE 5]. Nai THEN 

σετ ον στ ον ο ον ια بل ۷0 ال ول لد کلم‎ 


ENDIF 
CONTINUE 
CONTINUE 
WRITE (#, 205) %seeee R MATRIX ",KK,* X Ημ.’ HHHH 
DO 124° πα πα 
DO les LG 5 lN 
WRITE (€*, S000.  CRICIS و( ول‎ ο τα κα) 
CONTINUE 
WRITE (*,210) - 
CONTINUE 
WRITE (*, 205) °seeee SIL MAT RIX %%,KK,° X ',KH,? weww) 
Dig 126G IJ I HK 
ΘΕ; IL IM 
WRITE 6*5 300) τοις, τ. 
CONTINUE 
WRITE (+, 2109) 
CONTINUE 
WRITE (+, 205) ° wiae SZ m Rñ T R I X AK,’ X hie A469 
DO 128 1 = 1, KK 
98 1232 Ia Ia 
WRITE (FF, SOO ο πιο. Γι 
CONTINUE 
WRITE (#, 210) 
CONTINUE 


3x33 QU TPUT THE Z MATRIX «s 
WRITE (*,209) "sx Z MATRIX C,HKK,' X ^,KH," Qc»? 
WRITE (#, 212) 
DO 130 I-II Ue 
WRITE (X. 200) «σοι ον, Ὁ το hn 
WRITE (#, 210) 
CONTINUE 
WRITE (*, 2135) 


WRITE (G*, 419) 
READ  (*,200) ANSWER 
IF ((ANSWER .NE. 'Y') .AND. (ANSWER .NE. ’y’)) GOTO 21 


33339: ASK THE PARAMETERS FOR THE GRAPH 414 
WRITE (*, 210) 
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264 
eso 
EES 
267 
8ه‎ 
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ae 
27s 
274 
IB 
2/6 
er 
278 
3 
280 
281 
eo 
esas 
284 
کا‎ 
256 
287 
288 
2993 
0 
231 
eae 
eos 
294 
d 
236 
£937 
eB 
= 33 
200 
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202 
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3204 
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τό Microsoft FORTR 
WRITE (+, +*+) Tumeur EN TER PLOT FARAME 
WRITE (#, 410) 

READ (*, *) AZIM 

WRITE (*,411) 

READ (*, *) ει εν 

WRITE (#*, 413) 

READ (3, *) ITRIM 

WRITE (*, 414) 

READ (3, 3€) IDIV 

WRITE (*,41t5) 

READ (*, 199) CTEXT 

WRITE (*, 451) 

READ (*, 200) ANSWER 


AN7 


N 
ER 


Mi >i 
3) 
U 


3-33 INITIALIZES PLOT s 

IF ο ασ τα τα ο SER: ANSWER . EO. *y")) ΤΗΕΝ 
οαβι ρα σ ο ο 

= == ` 
CALL PLOTS (0, 39, 993) 

ENDIF 


CALL WINDOW((XLOL, YLOL, XUPR, YUPR) 


3333. DRAW. THE. MESH SURFACE ΟΕ THE GRAPH sx 
ΕΕ: MESS OL SE πώ ντ ον οσα εἰ Εν ο. SS O. SS. ES, E.S, DIV, O. 
em OROTZ EONTR MASE VERTEA) 


3333 ANNOTATION OF THE GRAPH sx 

ΕΠΕ SYMBGOL(CI.O, 6. 723, O0. 29, CTEXTS, O. 0$, 20) 

C SW S y IE: 6. ظ 16 د قو ف کت د له (ا بل نت وا‎ ο. ο, 19) 
ει 5٧9 τος ο ο ο ο ο AZIMUTH: nuUo) 
CALL NUMBER (939. 0, 9993. 0, 0. 23, AZIM, 0. 0, 2) 

εως ο πες (5 P DID 9.2 ΕΕ Πτι ο ο 16) 
SAHE NMOMSER CEEE OR SSS 0) Oo S9 AUV ο. ος 5) 

DY = στις 1) / SO ὃν * ELEV 

CAEL BEDDING lU ZCI) DV XR YR) 

GAEL ΥΜΕ πα ο کت‎ O OD) 

EEE SMEG شت لا اول 0 9 .لل‎ ο SUMIGIN',O90.0,10) 


eee OUTPUT THE GRAPH «κε 

CALE ο ο το د(‎ 39391 

WRITE (#, 416) 

READ (#, 200) ANSWER 

IF (CANSWER .E€Q. DS .GR. ANSWER CEO: ‘*y’)) GOTO هت‎ 


WRITE (*, 418) 
READ(*,200) ANSWER 
IF ((ANSWER .EQ. 'Y') .GR. (ANSWER .ED. “y’)) THEN 


tH ee! 


Heee FILL 9's THE TWO DIMENTIONAL GRID OF CONTROL POINTS τιν 


DG ise — = 1,26 


LSL 


fu fue ο 


(à fo [à fà [0 fü fu Q F 6b P F £6 £6 F t + O [u [D [u +" 


Lines 
207 
308 
203 
ت‎ 1 
الات‎ 
She 
E 3 
314 
ا‎ 
eve 
.لات‎ 
“18 
ο] 3 
GE 


RN 
اب‎ 
μαι 
<? k 
ed 
d irm td 
-— 
254 
saec 
لات‎ 
ed 


سا 


cu 


2 
2 
1 
7 MI EES E ΞΕΕ ΕΙ VII 
DD 12 go 
ZF(I,J) =O. OÓ 
11 CONTINUE 
ZFMAX = -9. 3E20 
ZFMIN = 9. 9E2090 
DK = (KK — 1) / (ا.ت‎ 
P= 5.141592 
DO 123 m = 1,KK 
DO 133 N = 1, KK 
RLPART = 9.9 
IMGPART = 0.0 
DO 134 L = 1,KK 
DO i134 K = 1,KK 
Ri¢1) = COS (-2*P+* (L-1) *(M—DK=-1) /KK) 
R1(2) = SIN(-2*P* (L=—1) * (M—-DK—-1) /KK) 
RzO1)0 = COS (-E*P* (K-11) *(N=—-DK=—1) /KK) 
R2(2) = SIN(-Z*P*(HM—1)*(N-DA-—-1)2)/HKH) r 
RLPART = RLPART + Z(L,K) #(R1 (1) *REC1L) 
* -=R1(Ə)*wRƏS (S) ) 
IMGPART = IMGPART + Z(L,K)*(R1(1)*RZ (2) 
* +*R1 (2) *¥ReEC(1L)) 
134 CONTINUE 
ZF (M,N) = SQRTC(RLPART##2 + IMGPART#*2) 
IF (ZF(M,N) .GT. ZFMAX) THEN 
ZFMAX = ZF(M,N) 
ENDIF 
IF (ZF(M, N) ΙΤ. ZEMIN) THEN 
ZFMIN = ZF(M, N) 
ENDIF 
Ds CONTINUE 
eeeee OUTPUT THE ZF MATRIX ***x«** 
WRITE (*,ZzO02) +++ FOURIER TRANSFORMATION °*,KK,’ X δν πμ. 
WRITE (#*, 212) 
DO 155 = 1,KK 
WRITS (*, 200) (Ze laa 
WRITE (#, 210) 
tsa CONTINUE 
WRITE (#5203) 
WRITE (*, 419) 
READ (*#, 200) ANSWER 
IF ((QNSWER .NE., *Y*') ΝΟ. (ANSWER ΝΕ " (('پ‎ GOTO 2 
*+*x** ASK THE PARAMETERS FOR THE GRAPH «κα 
aU WRITESCE, 210) 
WRITE (*,*) °*#* EN TER PLOT PRRAMETER S 
WRITE (*, 410) 
READ (3, 9€) AZIM 
WRITE (*, 411) 
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D Lines 
E 
So 
360 
261 
soc 
3e3 
364 
322 
266 
267 
368 
269 
370 
y E 
DN ds 
a7 
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DE 
276 
Da 
EXE 
279 
380 
381 
285 
285 
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1 386 
287 
388 
289 
Ego 
ga 
ms 
S32 
234 
i 
2285 
com 
298 
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400 
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402 
403 
404 
405 
406 
407 
408 
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a 
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03-35-95 
a mese 
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READ (*, *) E 
WRITE (#, 413) 

READ (4. *) ETRIM 
WRITE (#, 414) 

READ (3, 9) IDIV 
WRITE (#, 415) 

READ (*, 199) [IEXI 
WRITE (#, 451) 

READ (#, ZOO) ANSWER 


3*3 INITIALIZE PLOTEAG sere 

IF Ελ 5٣۴ EBS SYD ΟΝ. (ANSWER .EG@. yl) THEN 
CALL PLOTS (U, O, 5) 

ELSE 
BELL, PONS CO EE EE 

ENDIF 


WRITE (*, 4320) 
READ (*, 200) ANSWER 


CALL WINDOW(XLOL, YLOL, XUPR, YUPR) 


WE HANSER BRs VANYA Ra CANSHSR SEG. *y’)) THEN 
DLEV = (ZFMAX-ZFMIN) /FLOAT (KK) 
E AEE REE VEE THE SO O KET EES DEE ELE EES 
DO σου. => از‎ KKE 
LDIG(I) = 2 
LWGT(I) = 1 
CONTINUE 
S m Ὅἵ-ἵὨὉυ-υαποαοοσθ. 22 δ. Sy ZLLEV, LDIG, LWGT, 
IUS REIS, 10. 10) 
CeCe EOE tsea. 2070.2, CONTOUR ΜΑΡ" ο. ο 19) 
ELSE 
κκ κκ DRAW THE MESH SURFACE OF THE GRAPH κκ 


CALL MESHS (ZF, 26, 26, KK, KK, AZIM, ELEV, O. 5, 0. 5, 8. 25, 6. 5, IDIV, Ò, 


3, IPROJ, 1, ZLOW, 3, ITRIM, MASK, VERTEX) 


αλ ANNOTATION OF THE GRAPH sss»* 
CALL SYMBOL (5. 5, 2. 3, O. 2, AZIMUTH: O OOD 
CALL NUMEER (9399. 0, 399. 0, 0. £, AZIM, O. O, 9) 
ΓΕ ΞΕ ο ο ο ο ο ο ELE VATTON EMON O 10) 
ΕΓ: Naeem 29401), Jo. 0,0. 2c, ELEV,O. 0, 29 
DY = (ZF(1,1)/90.0) * ELEV 
EAL sDebD (tea. t., 2801,17 DV, XR, YR) 
ΕΓ. τος CCR MRS SEL 10:0, 19 
EAE SEO IO OAT OS ETT = ERIGIN',O.0,10) 
ENDIF 
ΓΕ. ΠΕΙ ο O επ ο Εν OO 29) 


CALL SYMBOL (6.0,6.5,0.2,°2-D DFT’, 0.9, 7) 


eeeee OUTPUT THE GRAPH +14 % 
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451 
455 
433 
424 
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436 
437 
438 
439 
440) 
441 
442 
44 3 
444 
445 
446 
447 
448 


Name 


ANSWER 
AZIM 
ΕΞΞ 

E EXT 
DK 
DLEV 
DY 
ELEY 


1:2 
zoo 
=O 
ipu 
eit 
ate 
cul ت‎ 
200 
4OC) 
451 
4.401 
402 
403 
404 
405 
406 
407 
408 
409 
410 
411 
413 
414 
4135 
416 
417 
418 
419 
429 


Tyne 


CHRR*1 
REAL 


CHAR* ZO 
REAL 
REAL 
REAL 
REAL 


7 Microsoft FORTROQN77 
CAGE αν ο ου 3) 
WRITE (*, 416) 
READ (*, 200) ANSWER 
IF ((ANSWER .EQ. *Y*')0) .QR. {ANSWER .EQ. *y’s) GOTO 
ENDIF 
WRITE (*, 417) 
READ (3*, 200) ANSWER 
IF ( (ANSWER .EQ. 'Y*) .QR. (ANSWER .EQ. *y’)) GOTO 10 
Siar 


FORMAT (AZO) 

FORMAT (A) 

FORMAT (/, 18X,AZ9, IZ, AS, Iz, AB, /) 

FORMAT () 

FORMAT (/, SX, 046) 

FORMAT (/, 2X, ' (AZIMUTH 320.0)", 46X,’ (AZIMUTH 230.0)!,/) 
FORMAT (/, 2X,’ (AZIMUTH O50.0)', 46X,’ (AZIMUTH 140.0)',7) 
FORMAT (10(F7. 2, i1X)) - 
FORMAT (9X, V) 

FORMAT (/, 5X, ' SEND GRAPH TO THE PRINTER(Y ar N): °,\) 
FORMAT(/,5X,’ NUMBER OF HORIZONTAL STATES(Nsi1to4): ',\) 
FORMAT (/,S5X, NUMBER OF VERTICAL STATES(Ms1t24): ',N) 
FORMAT (/, 5X, ' DIMENSION OF OUTPUT(itozE): *,\) 

FORMAT (SX, Al, 12, AS, 12, AZ, \) 

FORMAT (SX, AS, rz, AS. 12, AS, \) 

FORMAT (SX, AZ, 12,84, \) 

FORMAT (Sx, AS, 15. £2, AS, N) 

FORMAT (SX, AS, 12, AS. \) 

FORMAT (SX, AB, \) 

FORMAT (/, 5X, ’AZIMUTH(O.0 to 360.0 DEGREES): ',N) 
FORMAT (/,EX,' ELEVATION(30.0 t -90.0 DEGREES): °,\) 
FORMAT (/, SX,’ TRIM(O=NG, 1=%s, 2=Y¥s): ' V) 

FORMAT (/,5X,’2, 4 OR 8 SUBGRIDS: ',N) 
FORMAT, SX, TITLE OF GRAPH(UP TO $0 CHOR»): > N) 
FORMAT (/, 5X, "DO YOU WANT TO CHANGE PARAMETERS? ',^) 
FORMAT (/,5X,'DO YOU WANT TO REPEAT THE PROCESS? 1 
FORMAT (/,SX,’DO YOU WANT FOURIER TRANSFORMATION 7 ',N) 
FQRMAT(/,SX,'DO YOU WANT TO MAKE GRAPH ? *,\) 

FORMAT (/,5X,’DO YOU WANT CONTOUR MAP ? ',N) 

END 


Offset P Class 


32424 
22436 
INTRINSIC 
33448 
33488 


Sage 


33468 
33440 
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vse 


iQ 


Jane 9 
09-76-87 


13; 


κ. 
20 


See 


02/84 


ντι. 1 7 
FLOAT 

1 (INTEGER. 
IDIV INTEGERS 
EST INTEGER*Zz 
IMGPAR REAL 
IPROJ INTEGER 
IURI" Ἱπτεσςπακὸ 
IV REAL 

J 1۱١۱١ »نات‎ 
TT INTEGER *2 
K INTEGER*Z 
KK IN EGE REZ 
INTEGER*Z 
LDIG INTEGER*+2 
HJE INTEGERS 
LWGT INTEGERS 
m INTZGER*Z 
MASH INTEGER*Z 
N INTEGER te 
NRNG INTEGER*2 
ον REAL 

P REAL 

R REAL 

R1 REAL 

R2 REAL 
RLPART REAL 

ET REAL 

S2 REAL 

SIN 

SQRT 

TEMP REAL 
TEMP1 REAL 

TRM REAL 

REAL‏ لا 
VERTEX REAL‏ 

XEGE REAL 

XR REAL 

XUPR REAL 

"πω. REAL 

YR REAL 

YUPR REAL 

i REAL 

S REAL 
ZMAX REAL 
ZFMIN REAL 

το) REAL 

ZLOW REAL 


23158 
33446 
پس پس پس —- 

ت ت ت ت 

Fara ده‎ = 
-ω« δ 
στα 
334544 
33042 
33576 
323454 
سه‎ πον ړپ سم خم‎ 


. د 


cc 
So oe 
33134 
ت‎ 
22284 
06+ 
33164 
S616 
su 
23160 
33090 
33432 


320Z2& 
332034 
332308 
10818 


oc 
11616 
ulis 
33472 
33146 
33142 
33476 
alt 
O 
2704 
32480 
22484 
$408 


55154 


ΙΝΤΗΙΝΒΙΌΟ 
/WORM / 
/WORM / 
/WORM / 


THR INS iC 


IN RINSE 
/WORMK / 
/WORK  / 
WORK / 
/WORK / 
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Micrasort -ORTRONZ? 


D Line# 1 


Name Type 


MAIN 
ME SHS 
NUMBER 
PSDED 
PLOT 
PLOTS 
SYMBOL 
WINDOW 
WORK 
ZCNTUR 
ZLEVEE 


Pass One 


11680 


Class 


PROGRAM 
SUBROUTINE 
SUBROUTINE 
SUBROUTINE 
SUBROUTINE 
SUBROUTINE 
SUBROUTINE 
SUBROUTINE 
COMMON 
SUBROUT INE 
SUEROUT INE 


Na Errors Detected 


448 Saurce Lines 
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APPENDIX E 


ase 1 
Oder حدم‎ 
2 
سا‎ lias 1 7 Mrercnegrs 20347738 77 V5.0 U2/84 
1 “SLARGE 
2 $STORAGE: = 
SS SAGES LESS 
ج‎ 
SE Ji AE AE He E HE JE -ᾱ- -ᾱ- -ᾱ- -ᾱ- -ᾱ- -ᾱ- -ᾱ- -ᾱ- -ᾱ- هي خچ چي‎ -ᾱ- -ᾱ- چي د چ ې د‎ HE HE DIR EAR HE HE HB چ‎ HE de JE AE EE چب چب ې‎ οᾱ- -ᾱ- -ᾱ- -ᾱ- -ᾱ- -ᾱ- -ᾱ- -ᾱ- -ᾱ- -ᾱ- -ᾱ- -ᾱ- --- -ᾱ- ې خي پچ يد‎ ee 
5 c * * 
peu e THE PURPOSE OF THIS PROGRAM IS TO CODE TRE {=D (DISCRETE * 
cm + TIME) SYSTEM TO A Z2-D SPACIAL SYSTEM. * 
?2 C * * 
ΠΕ * EVANGELIS THEOFILOU + 
ήτο db db db d ecce Gh ie n de de de de ae ناک‎ -ᾱ- -ᾱ- P de am de Poe db de AB TO db db e M de eH dI db b e UR UI d A -ᾱ- -ᾱ- -ᾱ- -ᾱ- -ᾱ- -ᾱ- -ᾱ- -ᾱ- -ᾱ- -ᾱ- سي چب ما اک ښپ هي‎ -ᾱ- خا ني ا ا‎ 
13E PROGRAM 2D-DATA-=FIELD 
1.5 
14 C πες VARIABLE DECLARATIONS ** 
Ll REAL ο το ο  Ὁ σοι Sine) ys NL Cel Me Ce), PRM, GO), IVC SO). 
1 * IMGPSRT 
E CHARACTER#L ANSWER 
18 
T9 C Heee VARIABLI DECLARATIONS FQR PLOT38 «444+ τ 
= CHARACTER+=20 CTEXT 
=. COMMON ANRA EE E at S SAS T OR CECO EEE CES, 
2 * LDIG(26) , LWGT (26) , MASK (2000) , VERTEX (186) 
> 
z DATA Meee. O07, (Leo. O07, XUPR7S. S/, YURPR/ 7.0/7, 1 
25 * ZLOW/ 1. OESS/, IPROI/0O/, NRNG/100/ 
2 = و نه وه په‎ CHR FR ROG GR HH MAIN P RIGRAM Ad dE JE JE te He ها ها ې‎ ee 
ΞΕ eeeee OSK THE SEQUIRED VALUES FOR THE MODEL *e 
20 
ik 2 πίτες (* 305) 
oe READ (æ, +) M 
SS NEMESIS ο, ἂν OT CSI GOTO δ 
24 Ὁ WRITE (+, 404) 
= READ (*.*) M 
5 DNE ο ο a 661. wold) ΞΟΤῸ 5 
Sm 
du wRIGE (æ, 401) 
22 READ (#*, 290) ANSWER 
40 LF ο ο τω ۷٢ ( NM ANSWER τω, wv*)) THEN 
41 C «κος SILL O'S THE TWO DIMENTIONAL GRID CF CONTROL POINTS *«-- 
هه‎ π᾿, Oat eee o 
z 4 ت‎ DU OG NEN. 29 
2 44 ο ου τς ὁ 
2 = 36 CONTINUE 
46 
47 C ον ENTER VALUES FOR Y MATRIX se 
48 DO 37 I = 1, M*N 
1 43 ARITE(*,408) 'Y(C .وژ‎ 
1 eo EAD LC, FJ ROL, 1) 
i SI au CONTINUE 


1517 


14 


ace = 
5-5 
τες εν 
D wines 1 7 MIE SSDS FORTHEAN7Z?7 VE ES MEZA 
== ENDIF 
= M ((RMNSWER ασ Y) DR. (ANSWER EG DO GOTO 4 
54 
= WRITE (æ, 402) 
= READ (*, 200) ANSWER 
aia τα CCANSWER JEG, *Y°) ΠΠ. (ANSWER ١ ۰ (١٢ τα. 
EE νε INITIALIZE THE TRANSITION MAOTRIX sx 
= 2 BO 38 t = ifaw 
1 60 DO 38 J = 1, M+N 
= 61 TRM(I, J) Ξξῶ. ὁ 
z ez 38 CONTINUE 
63 
64 C 4*3 ENTER VALUES FOR TRANSITION MATRIX «ας 
= DO 93 I = 1,m+N 
1 s Dj 99.3 94 1 WM 
= 6 7 WRITE ORAO ΕΤΕ, ΑΗ 
= 63 READ (*, *) TRM(CI,J) 5 
= 22 22 CONTINUE 
70 END TS 
τα 
τοι ας 3*3 INITIALIZES R AND S ARRAYS s 
το DE ROG ERRES 
1 74 DO dO JE OS ° 
2 E με. ΞΕ εὐ 
2 yE QUIS uU 
2 YY 100 CONTINUE 
z 
vows e 3 INITIOLIZE INPUT VECTOR w***+- 
80 DO σι τι = ese 
1 9p ντο το. ο 
1 = 101 CONTINUE 
BS 
84 WRITE (#, 211) ‘ENTER INIT AL CONDITIONS ΕΓ ΠΙΕΙ’. 
Sa DG 102 1 = 1m 
1 86 σι τς πε ο X n ee 
1 87 READ (κκ...) RCI, 1) 
83 102 CONTINUE 
89 
ΠΩ WRIIE (*, 211) “ENTER INDE CONDI) LenS FOR ven ee SH’ 
En DO 193 I = 1,N 
1 im ΠΕ τ“, S ID : نت‎ 
EX READ (gero CI 


34 103 CONTINUE 


= 
- 


36 WRITE (*, 211) ENTER VALUES FOR THE INPUT VECTOR" 

37 IV(1) = 1.0 

38 WRITE (#*, 406) "ao1: ° 

393 READ (*,*) IV(M+1) 

100 

bo IF ((AMNSWER .ED. *Y*) .OR. (ANSWER .EQ. *'y*)9 GOTO Ss 
102 G κε INITIALIZE TRANSITION MATRIX HH ¥ 
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139 
140 
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147 
148 
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7 Micrasaft FOüWTRON77 V3. 


DO 194 I = 1,“ 
DO 104998 1,25 
TRM(I, J) = Q.Q 
104 CONTINUE 
WRITE (#*,211) ‘ENTER ZLEMENTS OF THE TRANSITION MATRIX? 
TRM (M+1, MN) = -IV(M+1) 
WRITE (Æ+, 406) Tall: ° 
READ  (*,*) TEMP 
TRM(1,M) = -TEMP 
τισ ο = -1.0 
WRITS CC DS :الد‎ ὉἹ 
READ εφ) TEMP 
TRM(M+1, M) = TEMP - TRM(1, M) * TRM(M+1, M+N) 


DO 105 = = ۷٨ 
TRM(I,I-1) = 1.0 
105 CONTINUE 


DO 106 I = 2+m,m+N 
TRM (I,I-1) = 1.6 
106 CONTINUE 


dc UNIS 
DO 107 I = 1, Næm : 
DO 108 J = 1, M+N 
te eee. ©) THEN 
DO 109 JJ = 1, M+N 


TP 2s. Ὁ πι τιν = RtJ.1+1) + 
* elo © Ths) 
Pecos GT. πμ πο σα τιν - 
* Sto I=", leer RM CS, J) 
109 CONTINUE 
πι τι = atl, ivi IVS) eu 
ENDIF 


IF (J .GT. Μ) THEN 
DO 110 JJ = 1,M+N 
IF (JJ .LE. M) S(J-M, I-1) 2 S(Q-M,I-1) - 


* ROJ, I) TRMCI, JJ) 
IF (JJ .GT. M) S(J-M,I-1) = S(J-M,I+1) + 
< S(JJ—-M, I) *TRM (J, JJ) 
110 CONTINUE 
S(J=m,I+1) = S(J=m,I+1) - ΙΝ (21) κ!) 
ENDIF 
108 CONTINUE 
U = 0.90 


107 GENTINUE 


WRITE (*,211) ?*s*«x«* INPUT VECTOR s»? 
WRITE CE SOO) (IVI), I = 1,M+N) 
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154 


co 
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1 D 
1 Ley 
158 

53 

160 

161 

1 162 
163 

164 

165 

166 

167 

1 168 
169 
170 
171 
172 
173 
174 
IX 
176 
177 
178 
1 د‎ 
180 
Pend 
182 
1 193 
1 184 
185 

186 

B 

188 

1 189 
1 190 
1 οι 
132 

135 
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